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ABSTRACT 

Various  methods  have  been  introduced  for  dynamic  braking 
of  ships.   Depending  on  the  prime  mover  and  transmission 
system,  dynamic  braking  is  accomplished  by  electric  power 
dissipation  (in  electric  driven  ships)  or  clutches  and  brakes, 
or  air  compression  (in  diesel  engine  propulsion  systems) . 

In  this  study  a  gas  turbine  -  electric  drive  combination 
is  used  as  the  means  of  propulsion;  the  usual  method  for 
dynamic  braking  is  the  electric  power  dissipation  method 
using  resistors. 

This  thesis  is  concerned  with  an  investigation  of  dif- 
ferent methods  for  dynamic  braking  in  an  attempt  to  find  a 
way  to  decrease  the  size  of  the  braking  resistors  needed  and 
their  associated  equipment,  or  to  eliminate  the  need  for  these 
by  using  another  type  of  dynamic  braking. 
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I.   INTRODUCTION 

This  thesis  is  concerned  with  the  study  of  a  ship's 
dynamic  braking  and  reversing  system  using  a  combination  of 
gas  turbine  -  electric  propulsion  motors. 

Because  of  the  ever  present  danger  of  collisions,  the 
problem  of  crash-stopping  a  ship  is  becoming  more  urgent  both 
in  the  Navy  and  the  Merchant  Marine. 

The  use  of  a  reduction  gear  mechanism  for  gas  turbine 
combatant  ships  is  quite  attractive  but  in  this  case,  during 
the  reversing  phase  it  introduces  a  complex  reduction  gear 
mechanism  requiring  large  clutches  for  the  accomplishment  of 
forward  and  astern  motion;  thus,  an  electric  drive  is  also 
attractive.   In  this  study  the  electric  drive  has  been  assumed 
and  various  methods  for  dynamic  braking  are  investigated. 

The  first  part  of  the  thesis  has  been  devoted  to  the 
study  of  different  open  loop  systems  in  order  to  select  the 
one  that  best  suits  the  required  specifications;  these  have 
been  defined  as  follows: 

a.  Stop  the  ship  in  minimum  time. 

b.  Try  to  minimize  additional  equipment  for  dynamic 
braking. 

The  second  part  of  the  thesis  is  concerned  with  the 
study  of  the  control  system  selected  in  the  first  part  imple- 
mented by  a  microprocessor-based  controller. 
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II.   THE  MODEL  SHIP 

For  the  simulation  of  the  problem  in  the  computer  (IBM 
360),  a  model  ship  was  derived.   For  this,  a  destroyer  escort 
type  vessel  was  selected  with  an  approximate  displacement  of 
4000  tons,  a  length  of  400  feet  and  a  top  speed  of  30  knots 
[Ref .  1] . 

A  single  fixed  propeller  was  selected  and  the  propulsion 
plant  consists  of  a  marine  gas  turbine  (FT4A-2) ,  which  directly 
drives  the  main  propulsion  generator  to  produce  the  necessary 
electric  power  for  the  main  motor  which  drives  the  propeller 
shaft.   A  schematic  diagram  of  the  propulsion  system  is  shown 
in  Figure  1. 

The  simulation  of  the  ship  and  propulsion  plant  uses  the 
above  model  with  somewhat  idealized  hydrodynamic  and  propulsion 
plant  dynamics. 

Since  the  equations  governing  the  response  of  the  ship 
and  propulsion  plant  are  complicated  non-linear  differential 
equations,  the  digital  computer  was  used  to  solve  the  equa- 
tions using  numerical  analysis  methods. 

As  a  consequence  of  the  idealized  model  the  following 
non-linear  functions  were  stored  in  the  computer  memory  as 
look  up  tables  during  the  simulation: 

a.  Thrust  reduction  coefficient       1  -  t 

b.  Wake  fraction  coefficient  1  -  W 

c.  Ship's  resistance  R 
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d.   Propeller  thrust  coefficient        C 
and   e.   Propeller  torque  coefficient        C 

q 

where  the  thrust  reduction  coefficient  (1-t)  is  a  measure  of 
the  propeller  and  hull  interaction,  and  the  wake  fraction  co- 
efficient (1-W)  is  a  measure  of  the  water  velocity  near  the 
propeller;  the  propeller  thrust  (CT) ,  and  torque  (C  )  coef- 
ficients  are  measures  of  the  efficiency  of  the  propeller. 

All  the  above  parameters,  along  with  the  corresponding 
computer  programs,  are  given  in  Appendix  B. 
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Fig . 1-  Schematic  representation  of  the  propulsion 
system 
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III.   DYNAMIC  BRAKING 

The  propeller  speed  can  be  reduced  quickly  from  some 
initial  large  valve  during  the  coast  down  phase  by  dynamic 
braking. 

Depending  on  the  prime  mover  and  the  transmission  system, 
dynamic  braking  is  accomplished  by  electric  power  dissipation 
(with  electric  drive) ,  clutches  (reversing  gear) ,  air  com- 
pression (diesel  engine) ,  or  with  a  water  or  air  brake  on  the 
shaft. 

During  dynamic  braking,  the  power  delivered  by  the  drive 
train  must  be  equal  to  the  power  absorbed  by  the  braking  system 
and  the  losses.   Since  the  drive  train  with  its  inertia  is 
being  decelerated  rapidly,  there  is  a  large  "transient  inertial 
horsepower"  due  to  the  shaft  and  propeller  speed  deceleration, 
as  well  as  rotational  power. 

In  terms  of  a  torque  equation,  this  can  be  expressed  as: 

2TT  Zt   ♦  ||E  =  Qmot  -  Qp  -  Qfr  -  QBR        (1) 

where,  Q    is  the  torque  developed  by  the  electric  motor  on 

the  propeller  shaft 
Q    is  the  torque  developed  on  the  propeller  by  the 

mass  of  water  passing  through  its  blades  due  to 

the  motion  of  the  ship. 
Q    is  the  friction  torque  on  the  shaft  and  motor 

FR 

Q,   ,    is  the  torque  developed  by  the  dynamic  braking 

device,  opposing  the  motion  of  the  propeller  shaft 
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I     is  the  total  moment  of  inertia  of  the  motor  rotor, 
propeller  shaft  and  propeller  (i.e.,  1=1   + 

shft     p 
N    is  the  speed  of  the  propeller. 

For  this  study  the  following  means  of  dynamic  braking 
will  be  considered: 

a.  Electric  power  dissipation 

b.  An  air  brake 

c.  Combinations  of  the  above 

d.  Use  of  the  inertia  of  moving  parts  as  load  and 
absorbing  medium. 

All  runs  have  assumed  an  initial  maximum  speed  of  30 
knots  and  a  maximum  propeller  speed  of  230  rpm. 
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IV.   SHIP'S  MOTION  EQUATIONS 

A  moving  ship  is  a  body  with  six  degrees  of  freedom. 
These  degrees  of  freedom  are  generally  chosen  as  follows: 

a.  Linear  displacement  along  the  three  axes  through  the 
center  of  gravity. 

b.  Rotations  around  the  three  axes  through  the  center 
of  gravity. 

If  we  assume  ship's  motion  in  calm  water  with  no  turning 
maneuvers,  then  the  motion  of  the  ship  is  best  described  by 
Newton's  law  of  motion  for  the  ship  in  translation  and  the 
propeller  in  rotation  [2],  as  follows: 

thrust  equation 

M~=(l-t).z  .T  (V  ,n)-R(V)      (2) 
at       p  p  p 

and   torque   equation 

24t.I.  4?p=  Q     +.  -  Q^     -Q  r^ 

t  dir      xmot       xfr     xp  (3) 

where       M  is  the  mass  of  the  ship 

V  is  the  velocity  of  the  ship 

Z   is  the  number  of  shafts 
P 

T  (V   )  is  the  propeller  thrust  as  a  function  of  the 
p  p,n  propeller  advance  speed  (V  )  and  propeller 
speed  (N) 

and         R  is  the  ship's  resistance  as  a  function  of  speed. 

The  input  data  required  to  simulate  an  acceleration  or 

deceleration  manuever  is  developed  from  the  characteristics 

of  the  ship  and  propulsion  system  being  modelled,  and  may  be 
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summarized  as  follows: 

Ship's  displacement 

Propeller  parameters 

Non-dimentional  propeller  performance  data 

Prime  mover  torque  characteristics 

Polar  moment  of  inertia  of  rotating  system 
The  coefficient  of  added  mass  (c) ,  representing  the  mass  of 
water  entrained  by  the  ship,  has  been  assumed  to  be  0.08  in 
all  cases  as  recommended  in  Ref.  12. 

Ship's  resistance  (R)  has  been  assumed  as  a  function  of 
ship's  speed.   In  addition,  the  hull-propeller  interaction 
is  also  assumed  to  be  a  function  of  ship  speed.   The  wake 
factor  (1-W)  and  the  thrust,  reduction  factor  (1-t)  are  gener- 
ally available  from  self-propelled  model  tests.   The  relative 
rotational  efficiency  is  assumed  to  be  constant  and  equal  to 
unity.   The  net  thrust  applied  to  the  ship  is  the  product  of 

the  propeller  thrust  (T  )  and  the  thrust  reduction  factor 

P 

(1-t) .   The  ship  speed  (V)  is  multiplied  by  the  wake  factor 
(1-W)  to  give  the  velocity  of  advance  of  the  propeller  (V  ) . 

In  terms  of  equations  the  above  are  given  as  follows 
[12]: 

Propeller  thrust  equation 

|  Tp=(l-t).Ta  (lf) 

where  TA  is  the  propeller  theoretical  thrust 

and 

Velocity  of  advance  of  propeller 
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V  =(l-w).V  (5) 

P 

where  V  is  the  velocity  of  the  ship. 

The  data  required  for  adequate  representation  of  fixed 
pitch  propeller  performance,  throughout  the  range  of  operating 
conditions  encountered  during  ship  acceleration  and  deceleration 
maneuvers  have  been  developed  by  Miniovitch  [2].   Miniovitch 
has  published  dimentionless  performance  data  for  three-bladed 
propellers  with  blade  area  ratios  from  0.6  to  1.6  in  incre- 
ments of  0.2. 

The  data  are  presented  in  the  form  of  thrust  (k  )  and 
torque  (k  )  coefficients  versus  the  advance  coefficient  (J) 

q 

2  2 

and  also  in  the  form  of  K  /J   and  K  /J   vs  1/J.   The  modified 

2  2 

thrust  (k  /J  )  and  torque  (k  /J  )  coefficients  versus  the 

reciprocal  of  the  advance  coefficient  (1/J)  are  required  in 

order  to  account  for  the  case  where  the  propeller  rotational 

speed  is  very  low  or  equal  to  zero.   The  data  extends  over 

the  four  distinct  operating  conditions  to  be  encountered  by 

the  propeller,  namely: 

Ship  travelling  ahead,  propeller  turning  ahead 

Ship  travelling  ahead,  propeller  turning  astern 

Ship  travelling  astern,  propeller  turning  astern 

Ship  travelling  astern,  propeller  turning  ahead. 

The  data  for  a  three  bladed  propeller  have  been  stored  in 

the  computer  as  a  look  up  table . 
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The  above,  in  terms  of  equations  can  be  expressed  as 
follows : 

Propeller  advance  coefficient 

J=N  .D/(V  2+(N  .D)2  )1/2  ,-. 

p    p    p  vb; 


Propeller  thrust 


T  =C„.p.D2.(V  2+(N  .D)2)  (7) 

a  T  r      p    p 


and  Propeller  torque 


Q  =C  .p.D3.(V  2+(N  .D)2)  (8) 

T  q  *         p   p 


where 

C_  is  the  propeller  thrust  coefficient 
C   is  the  propeller  torque  coefficient 

D  is  the  propeller  diameter 

p  is  the  sea  water  density 
N   is  the  propeller  speed 

The  polar  moment  of  inertia  of  the  rotating  system  has  been 
calculated  in  Appendix  A.   The  inertia  of  the  water  mass 
entrained  by  the  propeller  has  also  been  calculated. 

If  there  is  more  than  one  propulsion  shaft,  then  Eq.  (3) 
will  be  repeated  for  each  additional  shaft  [2] . 

The  direction  of  rotation  of  the  free  turbine  is  always 
clockwise  as  is  the  direction  of  rotation  of  the  generator 
set.   The  direction  of  rotation  of  the  motor  and  the  propeller 
shaft  is  assumed  positive  if  the  shaft  rotates  clockwise  when 
viewed  from  the  shaft  end. 
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Figure  2  shows  the  assumed  rotation  and  the  relation  of 
the  different  torques  acting  on  the  propeller  shaft  for  for- 
ward motion  of  the  ship. 


Fig. 2-  Assumed  notation  on  application  of  torques 
and  direction  of  rotation. 


The  generator  speed  (N    )  is  related  to  the  gas  turbine 


speed  (N   )  by  the  relation: 


N  =N 
gt  gen 


(9) 


The  frictional  torque  of  the  propeller  shaft  (QFR)  has  been 
assumed  constant  throughout  the  simulation;  actually  it  is  a 
function  of  the  shaft  speed  and  can  be  represented  also  as  a 
stored  table  look  up  [7] .   The  frictional  torques  of  the 
motor  and  generator  and  their  windage  torques  have  been 
assumed  negligible. 

A.    THE  GAS  TURBINE-GENERATOR  SYSTEM 

The  gas  turbine-generator  set  speed  is  controlled  by  the 
influence  of  the  engine  torque,  the  generator  load  torque  and 
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the  friction  losses  of  the  system.   This  can  be  expressed 

as  follows  [6] : 

dN 


2rt.I.    ,  .»t=  Qtt  -Q         -Q* 
tot  dt43         E       gen     xfr 


(10) 


where 


Q   is  the  gas  turbine  drive  torque.   It  is  a  non- 
El 

linear  function  of  fuel  consumption  (W  )  and 

r 

gas  turbine  speed  [8], 
Q     is  the  generator  load  torque  given  by  [4]: 


Q      =  K    .i-  (t).i  (t) 

xgen       g       fg  a 


(11) 


where 


where 


K  '  is  the  torque  constant  of  the  generator 

i,.  (t)  is  the  generator  field  current 
fg 

i  (t)  is  the  armature  current 
a 

I.  ,  is  the  sum  of  inertias  of  motor  rotor 
tot 

(I    )  and  gas  turbines  power  turbine  (I  . ) 

gen      3  r  gt 

Qf   is  the  friction  torque  which  is  made  up  of 
brush  friction  and  windage  components,  given 
by  [4]: 


xfr  xbr    w  gt 


(12) 


Q,   is  the  brush  torque 


K-j  is  the  windage  coefficient 

The  engine  is  set  to  run  as  a  nominally  constant  speed  engine 
with  a  proportional  governor.  Since  the  output  torque  of  the 
gas  turbine  is  a  non-linear  function  of  fuel  consumption  and 
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gas- turbine  speed,  data  for  the  output  torque  have  also  been 
stored  in  the  computer  as  a  look  up  table. 

B.    ELECTRICAL  TRANSMISSION  SYSTEM 

The  back  emf  generated  by  the  generator  and  the  motor 
are  given  by  [4] : 

E  =  K  .i-  (t).N  .  (13) 

g   g  fS     gt 

and       E  =K  .i-.  Ct).N  (14) 

m  m  tm    p 


where 

ij=    ,    i£  are  the  field  currents  of  motor  and  gen- 
fm   fg 

erator,  respectively 

Kg  is  the  generator  electrical  constant 

Km  is  the  motor  electrical  constant 

Em  is  the  back  emf  of  the  motor 

Eg  is  the  back  emf  of  the  generator 
Ignoring  the  effects  of  the  armature  inductances  of  motor  and 
generator,  for  maneuvering  conditions  the  armature  current 
will  be  given  by 

i  (t)=(E  (t)-E  (t))/R.  (15) 

a     a    m     t 


where 

R =  R  +  R  (16) 

t   m    s 


and  where 

Rm  is  the  motor  armature  resistance 

Rg  is  the  generator  armature  resistance, 
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V.   DYNAMIC  BRAKING  USING  ELECTRIC  POWER  DISSIPATION  METHODS 

A  braking  resistor  is  utilized  here  in  order  to  dissipate 
the  power  produced  by  the  motor  due  to  the  motion  of  the  shaft. 

When  the  CRASH-STOP  command  is  given,  the  main  generator 
is  disconnected  from  the  system  and  the  motor  is  left  to  rotate 
due  to  its  inertia  plus  the  applied  torque  on  the  propeller 
due  to  the  moving  mass  of  water  through  its  blades  due  to  the 
motion  of  the  ship. 

Thus,  the  motor  effectively  becomes  a  generator;  by 
leaving  its  terminals  open,  since  no  current  is  circulating, 
little  opposing  torque  will  be  produced.   If,  on  the  other 
hand,  a  resistor  is  connected  at  the  terminals  of  the  motor, 
a  current  will  circulate  and  thus  a  torque  opposing  the  motion 
of  the  motor  will  be  produced  and  the  propeller  speed  will 
drop. 

In  order  to  have  maximum  opposing  torque  on  the  propeller 
shaft,  the  maintenance  of  maximum  circulating  armature  current 
is  required;  but  since  the  voltage  generated  at  the  terminal 
of  the  motor  -  when  functioning  as  a  generator  -  is  variable, 
the  above  condition  can  only  be  met  by  a  variable  braking 
resistor  or  by  a  variable  motor  field  current. 

In  this  study  only  the  variable  resistor  method  has  been 
investigated.   Figure  3  shows  a  schematic  representation  of 
the  propulsion  system  after  the  CRASH-STOP  command  has  been 
initiated. 
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The  gas  turbine  can  be  held  at  load,  if  desired,  with 

the  simultaneous  addition  of  a  fixed  resistor  R  as  shown  in 

o 

Figure  3. 

The  equivalent  circuit  of  the  motor  is  shown  in  Figure 
4  after  t  =  0,  i.e.,  after  the  CRASH-STOP  command  has  been 
issued. 

Applying  Kirchoff's  law  around  the  loop  we  get 

V(t)=  E L(t)-i  (t).R.  -L  £|a(t)      (17) 
m     m    a    m   m  dt 

where 

Vm  is  the  voltage  generated  at  the  motor  terminals 

and  Lm  is  the  motor  armature  inductance 

Due  to  the  addition  of  the  braking  resistor  the  motor  now  is 

effectively  working  as  a  generator;  the  load  being  the  braking 

resistor  R_  and  the  prime  mover  being  the  torque  Q   produced 
a  p 

on  the  propeller  due  to  the  motion  of  sea  water  through  its 
blades . 

We  already  have  assumed  that  the  reactance  of  the  arma- 
ture circuit  is  negligible  so  that  Eq.  (17)  is  transformed  to 

V  =E  (t)-i  (t).R  (18) 

mm    a    m 

Applying  Newton's  law  for  the  motor-propeller  system,  we  get 
Eq.  (19),  i.e., 

2n.(I        +It    +1  .  _.)     ^Jp  =0        -  Q     -Q_  (19) 

mot     L       shft       dt^      Tnot     yp     yfr 


where  I      ,    =  moment  of   inertia  of  motor 

mot 

IT        =  moment  of   inertia  of  propeller 

Li 

I  h-p.!.  =  moment  of  inertia  of  propulsion  shaft 
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Fig.  3.   Schematic  representation  of  the 

propulsion  system  after  disconnection  (t   0) 
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Fig.  4.   Equivalent  circuit  representation  of 

the  separately  excited  DC  motor  at  t   0 
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and       Q   .  ,  the  torque  produced  by  the  motor  is  given  by 

0  *=  K»  .i.  (t).i  (t)  (20) 

where     K  '  is  the  torque  constant  of  the  motor, 
m 

As  can  be  seen  from  Eq.  (20),  since  the  armature  current 

has  opposite  direction,  Q   ,  will  be  negative  as  far  as  i_  (t) 

mot  rm 

remains  positive,  i.e.,  as  long  as  we  do  not  reverse  the  motor. 
Therefore,  the  right  hand  side  of  Eq.  (19)  is  negative  and 
the  motor  decelerates  faster  than  in  the  case  where  the  ter- 
minals were  kept  open  (i.e.,  Q    =  0) . 

From  Eq.  (20)  it  can  also  be  seen  that  the  bigger  Q 
is,  the  faster  the  deceleration  of  the  propeller  shaft  will 
be.   Thus  the  armature  current  must  be  kept  at  maximum  rated 
value  which  can  be  accomplished  by  varying  the  braking  resistor 

V 

In  order  to  calculate  the  rate  of  change  of  the  braking 
resistor  R^,  for  maximum  deceleration,  the  program  of  Appendix 
B  was  used  with  the  armature  current  held  constant  at  its 
maximum  rated  value.   The  values  of  the  braking  resistor  can 
then  be  calculated  from  Eq.  (18)  where  solving  for  R„  we  get 

RB__n^_a_m—  (21) 

a 

where 

V  =RD.I  (22) 

m  B  a 


The  calculated  values  of  the  braking  resistor  are  shown 
in  Table  I  and  have  been  plotted  versus  time  in  Figure  11. 
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As  can  be  seen  the  law  of  variation  of  the  braking  resistor 
is  not  linear,  therefore  an  approximation  was  made  by  fitting 
straight  lines  to  the  calculated  values. 

The  values  of  the  linearized  braking  resistor  are  given 
in  Table  I  also,  and  they  are  plotted  versus  time  in  Figure 
12. 

At  the  instant  the  shaft  speed  reaches  a  certain  value 
determined  by  the  maximum  loading  capabilities  of  the  system 
and  assuming  that  the  generator  field  current  has  been  reversed, 
the  system  is  reconnected  as  before  and  astern  maneuvers  are 
applied. 

The  armature  current  in  the  loop  will  be  given  by 

i  .(t)=(E  (t)-E  (t))/(R+R  )  (23) 

ast     g  m      m  g 

where  Eg(t)  is  the  back  emf  of  the  generator  given  by  Eq.  (13) 
and  Em(t)  is  the  back  emf  of  the  motor  given  by  Eq.  (14) .   Rm 
and  Rg  are  the  armature  resistances  of  motor  and  generator 
respectively. 

The  results  of  this  method  of  dynamic  braking  using  a 
linearized  braking  resistor  are  shown  in  Figure  13  through 
Figure  16,  along  with  results  that  would  have  been  obtained 
if  the  theoretical  variable  resistor  was  used. 

In  particular  Figure  13  shows  the  variation  of  ship's 
velocity  versus  time.   Due  to  the  fact  that  the  values  of  the 
linearized  braking  resistor  are  quite  close  with  the  values 


This  value  has  been  calculated  in  Appendix  A  to  be  13  rpm. 
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of  the  theoretical  one,  the  differences  of  speed  variation 
are  negligible. 

In  Figure  14  the  variation  of  propeller  speed  is  shown 
versus  time  for  both  cases.   The  abrupt  change  of  the  curve 
slope  during  dynamic  braking,  (for  0<r<1.4  sec.)  is  noticeable 
compared  with  the  change  of  slope  during  astern  manuevers . 

In  Figure  15  and  Figure  16  a  power  dissipation  diagram 
is  presented.   In  particular  in  Figure  15  the  power  dissipated 
in  the  braking  resistor  (PRB)  and  the  internal  resistance  of 
the  motor  (PRM)  have  been  plotted  against  time  for  the  case 
where  the  theoretical  variable  braking  resistor  is  used.   Also, 
on  the  same  graph,  the  percentage  of  dissipated  power  in  the 
braking  resistor  and  the  internal  resistance  have  been  plotted 
(PERC1  and  PERC2  respectively) . 

On  the  other  hand,  in  Figure  16  the  same  values  but  for 
the  case  of  the  linearized  braking  resistor  have  been  plotted. 
The  difference  in  the  percentage  of  dissipated  power  can  be 
noticed  since  now  a  linear  resistor  is  being  used.   Due  to 
linearization,  less  power  on  a  percentage  basis  is  dissipated 
by  the  braking  resistor.   Thus,  the  percentage  of  power  dis- 
sipated on  the  internal  resistance  of  the  motor  increases. 

Due  to  the  fact  that  with  a  linearized  variable  braking 
resistor  the  percentage  of  power  dissipation  in  the  internal 
resistance  of  the  motor  increases,  another  method  for  dynamic 
braking  using  a  fixed  braking  resistor  was  selected. 

In  this  case,  as  it  can  be  expected,  the  time  required 
to  stop  the  propeller  and  eventually  the  ship,,  increases . 
This  is  because  the  current  circulating  in  the  loop  is  no 
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longer  maximum;  therefore  the  opposing  torque  generated  by 
the  motor  (given  by  Eq.  20) ,  decreases  along  with  the  current. 

The  time  required  to  stop  the  propeller  with  a  variable 
resistor  was  found  to  be  1.8  seconds  while  with  a  fixed  re- 
sistor this  time  increases  to  28  seconds. 

The  results  of  the  method  where  we  use  a  fixed  braking 
resistor  are  shown  in  Figure  17  to  Figure  19.   In  particular 
Figure  17  shows  the  variation  of  ship  speed  versus  time  and 
it  can  be  noticed  that  the  ship's  speed  at  t  =  30  sec  is  still 
above  12  knots  while  from  Figure  14  the  time  to  stop  the  ship 
when  a  variable  resistor  is  used,  was  found  to  be  2  3.2  seconds. 

In  Figure  18  the  variation  of  propeller  speed  versus  time 
has  been  plotted.   The  largest  deceleration  rate  of  the  pro- 
peller speed  is  at  the  beginning  of  the  dynamic  braking  phase, 
since  then  the  circulating  current  is  still  maximum.   But  as 
time  advances,  the  circulating  current  drops,  therefore  the 
opposing  torque  produced  on  the  propeller  shaft  becomes  less, 
eventually  leaving  the  system  in  a  feathering  propeller  con- 
dition.  At  t  =  28  sec  the  slope  changes  due  to  the  application 
of  the  astern  motion  on  the  motor. 

In  Figure  19  as  in  Figures  15  and  16,  a  power  dissipation 
diagram  is  presented.   The  constant  percentage  dissipation 
curves  PERC1  and  PERC2  (for  the  braking  resistor  and  the 
internal  resistance  of  the  motor  respectively)  can  be  noticed. 
They  indicate  that  the  percentage  power  dissipated  in  the 
armature  resistor  is  constant  at  3%  while  the  remaining  97% 
is  being  dissipated  in  the  fixed  braking  resistor. 
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1.7 


1.72  10"2 

1.7136 

10"2 

1.5973 

lO"2 

1.3819 

lO"2 

1.1954 

lO"2 

1.0321 

lO"2 

8.8861 

lO"3 

7.6073 

lO"3 

6.4422 

lO"3 

5.3603 

lO"3 

4.3190 

lO"3 

3.3386 

lO'3 

2.4184 

lO"3 

1.5493 

lO"3 

0.0 

0.0 

0.0 

0.0 

Table  I.   Theoretical  and  Approximated  Values  of  the  Braking 

Resistor  RB  (ohms) 

TIME THEORETICAL APPROXIMATED 

0.0  1.72  10"2  1.7  io~2 

0.  I  L'llf,  \Ci~Z  1.6015  10"2 

0.2  ]  ..r,„7  !  m~J  1.5029  10~2 

0.3  1.3819  10"  1.4044  10~2 

0.4  1.1954  in~J  1.3059  10_2 

0-5  i.o?  -1   K.~J  1.2074  10"2 

0.6  8.8861  in"  ;  1.1088  10~2 

0.7  7.6073  in"  !  1.0103  10"2 

0.-  6.4422  in"  9.1176  10"3 


0-9  5.3603  10"J          8.1323  10_3 

1.0  4.31  qn  in"  '*  7.1470  10-3 

1.1  3.3386  in"'1  6.1617  10_3 
1-2  2.4  1 R4  If)"3          5.1764  10~3 

1-3  1  .rl4qi  in"3             4.1911  10"3 

1-4  O.n  3.2058  10"3 

1-5  0-0  2.2205  10"3 

1-6  0-0  1.5000  10~3 


0.0 
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The  CRASH-STOP  command  was  executed  in  2  3.2  seconds  when 
a  variable  linear  braking  resistor  was  used  and  in  t  »  35  sec 
when  a  fixed  braking  resistor  was  used. 
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VI.   DYNAMIC  BRAKING  USING  MECHANICAL  POWER 
DISSIPATION  METHODS 

Since  the  scope  of  this  thesis  is  to  try  to  minimize 
the  additional  braking  devices  required  for  dynamic  braking, 
it  may  seem  strange  to  consider  the  introduction  of  a  mechani- 
cal braking  device. 

But,  on  the  other  hand,  if  the  sizes  of  the  resistors 
found  by  the  application  of  the  previous  method  are  kept  in 
mind,  associated  with  the  problems  encountered  with  heat  dis- 
sipation then  it  will  be  found  that  the  use  of  a  mechanical 
brake  may  help  in  reducing  the  size  of  the  resistors. 

There  is  a  large  variety  of  mechanical  brakes  which  can 
be  used,  such  as  water  brakes,  friction  brakes,  or  eddy  current 
type  brakes . 

Clutches  and  brakes  are  used  to  control  the  coupling  be- 
tween two  shafts  by  means  of  relatively  low  power  signals. 
This  coupling  can  be  either  of  the  continuous  type,  where  the 
degree  of  coupling  varies  smoothly  with  the  magnitude  of  the 
control  signal,  or  of  the  on-off  type  where  the  coupling  is 
either  zero  or  maximum. 

A  brake  is  not  a  torque  converter  but  a  transmitter  of 
torque.   The  brake  controls  speed  by  permitting  slip  with 
respect  to  the  prime  mover  shaft.   To  accomplish  this,  the 
brake  must  absorb  energy  and  thereby  generate  heat  [3] . 

Heating  generally  establishes  the  maximum  power  level 
that  can  be  controlled  by  the  brake  by  limiting  the  maximum 
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slip.   An  analysis  of  the  heating  problem  must  include  an 
examination  of  the  duty  cycle  and  the  inertial  and/or  dissi- 
pative  characteristics  of  the  load.   Heat  dissipation  is  of 
maximum  concern  in  brakes  that  require  the  slipping  of  friction 
shoes  where  heat  is  localized  and  wear  is  high. 

Since  no  data  were  available  on  a  real  brake,  an  air- 
clutch  device  has  been  implemented  in  place  of  an  air  brake 
from  Ref.  9. 

The  amount  of  torque  that  a  brake  can  pass  without  slip- 
ping is  a  function  of  the  type  of  brake  and  the  net  air  pressure 
used  to  engage  the  brake. 

For  a  35  inch  diameter  brake  the  torque  is  given  by 

Q  =  4550. p  .     (lb-ft)  (25) 

where  P   ,  is  the  net  air  pressure  activating  the  brake  in 

net  r  3 

psi.   The  brake  glands  which  expand  the  brake  shoes  to  make 
contact  with  the  brake  face  expand  inwards  against  the  centrif- 
ugal force  of  the  rotating  axis.   Consequently  the  air  pressure, 
(P  ) ,  required  to  counteract  this  force  subtracts  from  the 
supply  pressure  (P  ) .   Thus  the  net  air  pressure  to  the  brake 

9 


glands  is  given  by 


P  .  =  P  -P  (26) 

net    s  c 


The  supply  pressure  typically  increases  upon  inflation 
at  the  rate  of  5  psi/sec  and  decreases  upon  deflation  at  the 
rate  of  30  psi/sec.   Thus,  for  a  typical  supply  pressure  of 
150  psi  and  with  t  in  seconds,  the  deflating  supply  pressure 
is  given  by  ' 
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P  =  150  -  30. t  (27) 

s 


and  the  inflation  supply  pressure  is 

P  =  5.t  (28) 

s 

The  pressure  required  to  counteract  the  centrifugal  force 
is  given  by 

P  =  5  +  (5.9  10"5  ).N  2  (29) 

c  p 

where  N  =  60  .  n  (30) 

P       P 


therefore  P  =  5  +  (3.54  10  3)  .n  2  (31) 

c  p 

Combining  the  above  equations  yields  the  following  expressions 
for  the  torque  supplied  by  the  brake: 
During  engaging 

Q„T  =4550  (5t-5-3.54  10~3  .n  2)  (32) 

XCL  p 


and  during  disengaging 

QCL  =  4550  (145-30t-3.54  10_3.n  2)        (33) 

Then,  from  Newton's  law,  the  torque  equation  for  the  motor- 
propeller  system  is  given  by: 


M  I  *  -%  -V  <a.  m) 
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Figures  20  and  21  give  the  results  of  using  the  above 
mechanical  brake  as  a  dynamic  braking  device. 

In  particular,  Figure  20  shows  the  variation  of  ship's 
speed  versus  time  and  Figure  21  shows  the  variation  of  propel- 
ler speed  versus  time.   Notice  the  small  change  in  the  slope 
of  the  curve  during  deceleration.   This  occurs  due  to  the  fact 
that  there  is  a  delay  in  the  application  of  the  braking  device 
and  also  that  the  torque  produced  by  the  moving  water  through 
the  blades  of  the  peopeller  is  still  quite  large.   At  t  Z   10 
sec  the  change  of  slope  is  due  to  the  application  of  astern 
torque  on  the  propeller  shaft. 

Since  the  use  of  a  mechanical  brake  alone  as  a  dynamic 
braking  device  will  be  associated  with  maintenance  and  heat 
dissipation  problems,  a  combination  of  the  above  mechanical 
brake  and  the  electric  power  dissipation  methods  examined 
earlier  will  be  attempted. 

This  could  result  in  the  smoothing  of  some  of  the  problems 
usually  encountered  with  the  stand  alone  application  of  each 
system.   Reduction  of  the  size  of  the  braking  resistors  or 
reduction  of  the  heat  dissipation  on  the  mechanical  brake  are 
anticipated  benefits. 

For  this  case  Eq.  (34)  is  modified  by  the  addition  of 
Eq.  (20)  resulting  in 

^n   =0       .    -0      -0.      -0_  (35) 


where 


^tot  at?  =Qmot  '%  -Qfr  "QCL 


Itot~Imot  +  shft      p 
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For  the  case  where  the  mechanical  brake  is  used  in  com- 
bination with  the  linear  variable  resistor  no  significant 
advantage  is  gained  as  is  shown  in  Figures  22,  23,  and  24,  and 
the  reduction  of  the  size  of  the  braking  resistor  is  insig- 
nificant since  most  of  the  dynamic  braking  is  accomplished 
by  the  variable  braking  resistor. 

In  particular,  Figure  22  shows  the  variation  of  propeller 
speed  versus  time.   It  can  be  seen  that  for  0  <  t  <   2    sec  the 
deceleration  of  the  propeller  is  primarily  due  to  the  variable 
resistor  since  from  Eq.  (32)  it  can  be  seen  that  there  exists 
a  time  delay  in  the  application  of  the  mechanical  brake.   This 
deceleration  of  the  propeller  is  almost  the  same  as  that  shown 
in  Figure  14  where  only  a  linear  variable  resistor  is  used 
for  dynamic  braking. 

Figure  24  shows  the  torque  developed  by  the  mechanical 
brake  during  dynamic  braking.   As  can  be  seen  it  is  not  sig- 
nificant as  to  cause  any  change  in  the  dynamics  of  the  system. 
Figure  2  3  shows  the  variation  of  ship's  speed  versus  time. 

On  the  other  hand,  the  combination  of  the  mechanical 
brake  and  the  fixed  braking  resistor  offers  significnat  ad- 
vantages as  far  as  time  requirements  are  concerned.   Thus  the 
time  to  stop  the  propeller  reduces  from  28  seconds,  as  in  the 
case  where  a  fixed  braking  resistor  was  used  as  the  only  means 
of  dynamic  braking,  to  6.8  seconds  with  the  addition  of  the 
mechanical  brake.   Similarly,  the  time  needed  to  stop  the  ship 
decreases  from  t  »  35  sec  to  approximately  25  seconds. 

Figure  27  shows  the  torque  developed  by  the  air  brake. 
It  can  be  noticed  that  for  the  initialization  of  braking 
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action  approximately  1.2  seconds  are  required.   Also,  when 
disengaging  should  occur,  it  does  not  occur  immediately  but 
some  delay  is  encountered. 
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VII.   DYNAMIC  BRAKING  USING  THE  INERTIA  OF 
MOVING  PARTS  AS  LOAD 

Due  to  problems  of  space  requirements  and  heat  dissipation 
when  electric  or  mechanical  power  dissipation  devices  are  used, 
an  attempt  was  made  to  find  a  way  to  minimize  or  eliminate 
these  problems. 

Considering  first  the  effect  of  ship's  motion  on  the 
propulsion  system,  it  is  relevant  to  note  the  various  mechan- 
isms present  in  the  electric  drive  during  astern  maneuvers. 

The  first  stage  in  going  CRASH-STOP  is  with  positive 
voltage  at  the  generator  and  positive  armature  current,  the 
current  being  reduced  at  a  fast  rate.   Since  we  want  to  apply 
astern  conditions  as  fast  as  possible,  the  way  to  stop  the 
propeller  is  to  eventually  reverse  the  armature  current. 

The  second  stage  begins  the  moment  that  the  armature 
current  reaches  a  negative  value;  that  is,  from  this  moment 
on  the  generator  starts  to  motor  from  the  propeller's  kinetic 
energy 

During  this  second  stage  Eq.  (10)  becomes 

«VV  ^;  qe  V  ^  C36' 

where  the  relation  for  Q    is  given  by  Eq.  (11)  and  for  Qf 
by  Eq.  (12) . 

As  can  be  seen  from  the  above  equation  the  gas  turbine 
will  start  to  accelerate  unless  a  way  is  found  to  reduce  the 
effect  of  the  reversed  operation  of  the  propulsion  system. 
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Fig.i:a  -  Diagrammatic  representation  of  the  propulsion  system 
dyring  CRASH-  STOP  manouvres  .  (a) -First  stage  , 
(b) -second  stage  , (c) -third  sta^e  . 


43 


One  way  would  be,  if  possible,  to  disconnect  the  gas 
turbine  from  the  generator  set  so  that  overspeeding  of  the 
generator  would  not  have  any  effect  on  the  gas  turbine.   But 
in  the  model  system  we  have  selected,  the  gas  turbine  is 
rigidly  connected  to  the  generator  set,  so  this  cannot  be 
done. 

Because  of  the  limitations  on  the  overspeeding  of  the 
gas  turbine  given  in  Ref.  5  other  ways  have  to  be  found. 
The  limitations  are: 

Power  Turbine  Speed 
Maximum  continuous  speed  3600  rpm 

Topping  governor  reference  speed        3744  rpm 
Over speed  trip  setting  3960  rpm 

Since  maximum  opposing  torque  is  needed  on  the  propeller 
shaft,  this  means  that  maximum  current  will  flow  in  the 
armature  circuit.   This  can  only  be  accomplished  by  an  effec- 
tive field  generator  current  control  according  to  the  relation 

E  =  E  -  33  (1)    Volts  (37) 

g        m 


But  since 


E  =  K  .  i„  (t)  .N  .  ,0Qs 

g    g    fg      gt  (38) 


substituting  Eq.  (38)  into  Eq.  (37)  and  solving  for  If   gives 

I*  (t)  =-Sr  ~33 


fg     N"1.K  (39) 

gt  g 


This  has  been  calculated  in  Appendix  A. 
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Since  the  armature  current  (I  )  will  be  maintained  at  its 

3 

maximum  rated  value,  it  can  be  seen  from  Eq.  (36)  that  the 
only  contribution  that  can  reduce  the  gas  turbine  overspeeding 
effects  is  the  term  Q_,  i.e.,    the  gas  turbine  output  torque. 

Thus,  an  attempt  will  be  made  during  the  second  stage  to 
reduce  the  gas  turbine  output  torque.   Since  Q   is  a  function 
of  the  fuel  flow  rate  (Wf)  and  the  gas  turbine  speed  (N   ) 
the  reduction  of  the  fuel  flow  rate  in  the  gas  turbine  was 
selected  as  the  means  of  doing  this. 

At  the  CRASH-STOP  command,  the  fuel  flow  rate  is  reduced 
to  that  corresponding  to  the  idle  speed  of  the  gas  turbine 
and  the  armature  current  is  controlled  and  held  at  its  max- 
imum value  by  the  field  generator  current,  according  to  Eq. 
(33). 

Then,  according  to  Newton's  Law,  the  torque  equation  on 
the  propeller  shaft  becomes 

This  permits  rapid  deceleration  of  the  propeller.   It  must 
be  noted  here  that  in  a  real  system  the  overspeeding  effects 
of  the  second  stage  would  be  additionally  lessened  by  the 
brush  friction  and  generator  windage  and  friction  torques 
which  at  this  stage  have  been  assumed  negligible. 

The  third  stage  begins  as  soon  as  the  generator  voltage 
is  reversed.   From  Eq.  (38)  this  will  happed  when  the  propel- 
ler speed  is  at  approximately  13  rpm;  then  astern  motion  of 
the  shaft  can  be  effectively  applied.   This  is  accomplished 
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by  reversing  the  generator  field  current  gradually,  according 
to  Eq.  (39) ,  so  that  no  overloading  of  the  armature  circuit 
will  occur. 

Concurrently  at  this  stage,  the  fuel  flow  rate  at  the 
gas  turbine  is  increased  gradually  to  cause  the  turbine  speed 
to  assume  its  normal  operating  value  (3600  rpm)  . 

When  ship's  speed  approaches  zero,  a  gradual  decrease 
of  the  generator  field  current  will  bring  the  propeller  shaft 
to  rest.   This  can  be  accomplished  either  manually  by  the 
operator  or  automatically  by  the  control  system. 

The  results  of  this  method  are  shown  in  Figures  28,  29, 
30  and  31.   In  particular,  Figure  28  shows  the  variation  of 
propeller  speed  during  the  CRASH-STOP  maneuver.   The  time 
to  stop  the  propeller  was  found  to  be  1.5  seconds,  i.e.,  in 
approximately  6  revolutions  of  the  propeller  shaft.   In 
Figure  29  a  plot  of  the  law  of  variation  of  the  field  current 
is  shown,  according  to  Eq.  (39)  in  order  to  have  maximum 
armature  current  in  the  loop  resulting  in  maximum  deceleration 
conditions.   Figure  30  shows  the  speed  of  the  gas  turbine 
during  dynamic  braking,  where  the  overspeeding  effects  of 
this  method  can  be  noticed. 

The  gas  turbine  speed  increases  up  to  3660  rpm  during 

2 
the  second  phase.    Then  phase  three  is  initiated  and  astern 

conditions  are  applied. 


2 

In  accordance  with  overspeeding  limits  given  in  Ref.  5, 
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In  Figure  30  it  can  be  noticed  that  due  to  the  fact  that 
the  fuel  flow  rate  in  the  gas  turbine  cannot  be  increased 
immediately  to  the  value  corresponding  to  normal  operating 
conditions,  deceleration  of  the  gas  turbine  speed  takes  place 
until  the  flow  rate  of  fuel  builds  up  to  the  required  amount, 
where  the  gas  turbine  speed  is  brought  back  to  normal  operating 
conditions . 

Finally,  Figure  31  shows  the  variation  of  ship  velocity 
versus  time.   The  time  required  to  stop  the  ship,  from  a 
maximum  speed  of  30  knots,  was  found  to  be  2  5.2  seconds, 
which  is  quite  compatible  with  the  case  examined  earlier 
where  a  variable  linear  resistor  was  used  as  a  braking  device 
as  is  shown  in  Figure  13. 
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VIII.   DISCUSSION  OF  RESULTS 

The  results  to  be  discussed  are  shown  in  Table  II. 

The  use  of  a  variable  resistor  as  a  braking  device  offers 
the  advantage  of  fast  deceleration  of  the  propeller  shaft  in 
1.7  seconds  from  an  initial  speed  of  230  rpm  to  13  rpm.   This 
amounts  to  stopping  the  propeller  shaft  in  approximately  six 
revolutions. 

Along  with  this  advantage  are  disadvantages,  however. 
First,  the  use  of  a  variable  resistor  can  only  be  associated 
with  a  combination  of  series  or  parallel  resistors  along  with 
a  voltage  controlled  timing  mechanism.   Large  circuit  breakers 
will  also  have  to  be  incorporated.   Space  limitations  and  heat 
dissipation  problems  are  additional  disadvantages  which  will 
add  to  the  complexity  of  the  system. 

On  the  other  hand,  the  use  of  a  fixed  resistor  eases  the 
problems  of  complexity  and  high  maintenance  costs  associated 
with  the  timing  mechanism  and  the  circuit  breakers,  but  still 
has  the  additional  space  limitations  and  heat  dissipation 
problems  though  to  a  lesser  degree  due  to  the  reduced  size  of 
the  resistor  bank  (fixed) .  But  the  decrease  in  complexity  is 
paid  in  an  increase  of  time  results;  that  is  the  time  to  stop 
the  propeller  is  increased  to  27  seconds  and  the  time  to  stop 
the  ship  from  maximum  speed  is  much  greater  than  35  seconds. 

The  introduction  of  a  mechanical  brake  gives  almost  the 

same  timing  as  the  variable  resistor  braking  device  does  but 

still  space  requirements  and  heat  dissipation  problems  are 

present. 
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As  is  shown  in  Table  II,  the  combination  of  the  variable 
braking  resistor  and  the  mechanical  brake  does  not  offer  any 
significant  advantage  in  time  requirements  and  it  almost 
doubles  the  space  needed  for  accommodation  of  the  extra  mechan- 
ical brake.   Similarly  the  reduction  of  the  size  of  the  resistor 
is  insignificant. 

On  the  other  hand,  there  are  advantages  to  using  a  fixed 
resistor  combined  with  the  mechanical  brake.   There  the  time 
to  stop  the  propeller  is  found  to  be  6.8  seconds  instead  of 
28  seconds  when  the  fixed  resistor  is  used  alone  or  9.9  sec- 
onds when  the  mechanical  brake  is  used  alone.   The  time  re- 
quired to  stop  the  ship  from  maximum  ahead  speed  is  found  to 
be  approximately  25  seconds  instead  of  approximately  50  sec- 
onds when  the  fixed  resistor  is  used  alone,  or  25  seconds 
when  the  air  brake  is  used  alone.   But  still,  in  this  case, 
the  disadvantages  resulting  from  the  use  of  heat  dissipating 
elements  still  exist  along  with  increased  maintenance  problems 
due  to  the  addition  of  the  mechanical  brake. 

For  the  last  case  of  dynamic  braking  examined,  i.e., 
using  the  inertia  of  moving  parts  of  the  power  turbine  and 
the  generator  rotor,  some  assumptions  were  made. 

First,  it  was  assumed  that  the  fuel  flow  rate  in  the  gas 
turbine  could  be  reduced  from  the  operating  condition,  to  that 
corresponding  to  the  idle  speed  of  the  gas  turbine,  instantan- 
eously.  Second,  complete  controllability  of  the  generator 
field  current  was  assumed,  so  that  constant  maximum  reverse 
armature  current  can  be  maintained. 
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Under  the  above  assumptions,  the  simulation  of  the  system 
gave  the  following  results: 

Time  required  to  stop  the  propeller  from  maximum  speed 
to  zero  rpm  equals  1.5  seconds.   Time  required  to  stop  the 
ship  from  30  knots  to  zero  equals  2  3.2  seconds. 

For  the  operation  of  the  system  no  extra  equipment  is 
needed  such  as  braking  resistors  or  brakes,  nor  any  heat  dis- 
sipating elements  in  excess  of  those  already  existing.   The 
cost  of  the  installation  is  not  increased  due  to  addition  of 
extra  braking  devices  and  the  cost  of  maintenance  remains 
almost  the  same . 

In  order  to  be  able  to  compare  the  systems  discussed, 
as  far  as  performance  is  concerned,  another  assumption  which 
was  made  for  all  cases  examined  has  to  be  taken  into  account; 
that  is,  it  was  assumed  that  astern  conditions  could  be 
applied  immediately  at  maximum  operating  conditions.   The 
effect  of  this  assumption  is  that  the  time  required  to  stop 
the  ship  will  be  increased  by  an  equal  amount  for  all  cases 
examined. 

Based  on  the  above  discussion  and  having  in  mind  that 
the  scope  of  this  thesis  was  to  investigate  the  existing  types 
of  dynamic  braking  and  specifically  try  to  minimize  and  possibly 
eliminate  the  need  for  dynamic  braking  resistors,  the  last 
system  was  selected  for  further  investigation. 

A  block  diagram  of  the  open  loop  system  in  shown  in 
Figure  5 . 
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IX.   CONTROL  SYSTEM  DESCRIPTION 

A.    GENERAL 

The  system  simulation  was  set  up  for  the  individual  sub- 
systems outlined  below: 

1.  Gas-turbine-generator  system 

2.  Electric  transmission  system 

3.  Generator  field  current  control  system 

4.  Propeller-hull  system 

5 .  DC  motor . 

Subsystems  1,  2,  4,  and  5  have  already  been  discussed  in  the 
previous  sections  and  there  is  no  need  to  describe  them  here 
again.  For  references  purposes  the  equations  governing  each 
system  will  be  repeated  here. 

Gas  turbine-generator  system 

The  gas  turbine-generator  set  speed  is  controlled  by  the 
influence  of  the  engine  torque  and  the  generator  load  torque 


and  frxction,  or 


where 


2nl  f£H  -Qgen-Qfr 


I     =1-+.     +     I 

gt        gen 


where  the  gas  turbine  drive  torque  is  a  function  of  the  fuel 
flow  rate  and  the  speed,  or 

QE  =  f  (wf  '  V 
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The  generator  torque  is  given  by 

Q    =K   .  I_  .  I 

gen    g    fg    a 

and  the  friction  torque  is  made  up  of  brush  friction  and 
windage  components  assumed  to  be  of  the  form 

Q.  =  Q.   =  K  .  N  2 
xfr    br    w   g 

The  engine  is  set  to  run  at  constant  speed. 
Electrical  transmission  system 
The  transmission  equations  are 

E  =  K   .1-   .  N. 
g    g   fg    gt 


E  =K  .  I-  .  N 

m   m    fm    p 


Ignoring  the  effects  of  armature  inductance  for  maneuvering 
conditions,  the  armature  current  is  given  by 

where  Rt  =   Rg  +  Rm 


Propeller  and  hull  system  and  motor 

The  propeller  dynamics  are  described  by  the  following 
equations : 

2n.Ip  |E  .  ^  -Qp  .  Qfr 
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where  I  =  I  + 1  *«. + 1 

p         m         shft         p 


,. 


and  o     =  K       .1-     .   I 

Til        m         fm         a 


The  load  torque  on  the  propeller  is  defined  by 

Q  =  C   .  p  .  D3  .A 

where  D  is  the  propeller  diameter,  p  is  the  local  sea  water 

density  and 

A  =  (V  L   +  N  .Dz  ) 
a     p 

The  coefficient  C   is  defined  by  graphical  data  against 
the  modified  advance  coefficient  S  where 

S  =  V  /  A1/2 

cl 

The  ship  hull  dynamic  effects  are  given  by 

M^-  =  T  -T 
dt     p    v 

where 

T  =  (^  .  p  .D2  .A 
P    l 

and 


°T 


f(S) 


and  T   is  the  ship's  resistance  as  a  function  of  the  ship's 
velocity. 
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B.  GENERATOR  FIELD  CONTROL  SYSTEM 

The   generator   field  coil  equation   is   given  by 

L      dlfg  =  V  (t)     x       .R 

f  g  dt  c  f  g       f  g 

where     L,.   is  the  field  circuit  inductance 

fg 

V  (t)  is  the  voltage  of  the  field  circuit 
c 

R-   is  the  field  circuit  resistance. 

fg 

A  closed  loop  system  to  control  shaft  speed  using  the 
generator  field  current  to  achieve  speed  changes  would  be 
appropriate  for  the  following  reasons: 

1.  Speed  control  with  a  percentage  accuracy  better  than 
the  percentage  system  losses  is  not  necessary  and  the  operator 
could  compensate  for  the  losses  if  required. 

2.  During  turns  using  rudder  control  it  is  not  desirable 
to  compensate  for  the  resulting  shaft  speed  changes. 

3.  The  prime  mover  has  been  designed  to  run  as  a 
constant  speed  machine. 

The  control  system  then  is  basically  a  field  current 
controller. 

C.  DEVELOPMENT  OF  THE  FIELD  CURRENT  CONTROLLER 

From  the  discussion  of  the  open  loop  system,  it  has  been 
observed  that  the  armature  current  must  be  constant  in  order 
to  apply  maximum  opposing  torque  on  the  shaft.   From  the  re- 
lation 

E  =  E  -  33  Volts 
g    m 


substituting  the  equation  for  E   and  solving  for  I_   gives 
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I.  =  (E  -  33)/K  .N. 
fg     m       g  gt 


Due  to  advances  in  solid  state  electronics  and  micro- 
processor chips  the  use  of  a  microprocessor  based  controller 
was  assumed. 

Comparing  a  microprocessor  based  controller  with  the  up 
to  now  commonly  encountered  field  current  controllers  incor- 
porating thyristor  devices  [7] ,  we  obtain  the  following: 

1 .  Control  Room  Space  Considerations 

Because  of  the  high  density  packaging  used  in  micro- 
processor based  controllers,  large  amounts  of  computational 
and  control  components  can  be  placed  in  a  relatively  small 
enclosure  to  reduce  control  room  space  requirements. 

2 .  Component  Count 

Microprocessor  controllers  use  "ti:r*S?  shared"  com- 
ponents and  thus  reduce  component  count.   The  same  circuit 
that  compute  a  3-mode  algorithm  is  also  used  to  output  out- 
of-limit  alarms  and  extract  square  roots.   One  microprocessor 
controller,  which  can  be  used  in  lieu  of  up  to  six  analog 
controllers  and  numerous  computing  relays,  operates  at  least 
four  final  control  elements. 

3.  Accuracy 

There  is  no  question  that  digital  devices  are  more 
accurate  than  analog  devices.   Using  a  15-bit  word  length 
microprocessor  uncertainty  in  any  addition,  multiplication, 
division  or  square  root  is  about  0.003%.   Digital  devices  do 
not  suffer  from  drift  problems.   And  with  proper  design,  the 
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noise  immunity  of  digital  devices  can  be  made  superior  to  that 
of  analog  devices. 

4 .  Reliability 

The  reliability  of  a  microprocessor  based  digital 
controller  can  be  made  much  higher  than  the  conventional  control 
equipment  that  it  replaces.   High  temperatures  are  avoided  in 
microprocessor  controller  designs  by  virtue  of  the  low  power 
consumption  of  solid  state  electronics.   Although  the  reliability 
of  a  single  microprocessor  controller  may  or  may  not  be  super- 
ior to  a  single  analog  controller,  there  is  no  question  that 
a  single  microprocessor  controller  is  more  reliable  than  a 
system  consisting  of,  let  us  say,  six  analog  controllers,  30 
computing  relays,  and  a  panel  full  of  logic  circuitry.   Use 
of  a  microprocessor  also  allows  a  number  of  "tricks"  to  be 
played  by  the  programmer  to  enhance  basic  reliability  even 
further. 

Self-test  features  can  easily  be  incorporated  into 
the  program  and  failure  of  a  self-test  can  be  used  to  alert 
the  operator  that  a  replacement  should  occur. 

5.  Flexibility 

With  analog  and  discrete  components  in  a  control 
system  a  change  in  control  strategy  normally  requires  a  change 
in  control  equipment;  these  changes  can  be  easily  implemented 
by  software  changes  in  a  microprocessor  based  controller. 

6.  Cost 

Generally  speaking  for  a  single  application  the 
cost  of  analog  devices  may  be  smaller;  but  for  a  number  of 
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applications  the  pertinent  cost  of  the  microprocessor  based 
controller  decreases  rapidly. 

A  simplified  diagram  of  a  microprocessor  based  controlled 
is  shown  in  Figure  6 . 


1L=*> 


Figure  6.   A  real  time  u-controller  system 
A  dedicated  microcontroller  receives  information  from 
the  input  sensors  and  it  is  connected  to  the  plant  so  that 
it  may  exercise  direct  control. 

Input  to  a  microcontroller  can  be  from: 
push  buttons 
limit  switches 
voltage  sensors 
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temperature  sensors 
pressure  sensors,  etc. 
The  output  of  a  microcontroller  can  be  directed  to: 
starters 
solenoid  valves 
stepping  motors 
displays 

servomotors,  etc. 
Figure  7  shows  a  typical  programmable  controller  processor 
diagram  along  with  the  associated  units. 


( 

Figure  7.   The  programmable  controller/processor  system 

The  sensing  devices  of  the  process  are  connected  through 
converters  to  the  microcontroller.   These  converters  (trans- 
ducers) reduce  the  size  of  the  actual  signal  and  then  the 
signal  is  passed  through  A/D  converters  which  change  the 
analog  signals  into  digital.   After  processing,  the  output 
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data  are  passed  through  a  D/A  converter  which  transforms  them 
into  analog  signals  to  be  used  as  control  commands. 

In  the  actual  system  to  be  used  the  following  inputs  are 
required: 

a.  gas  turbine-generator  speed 

b.  fuel  flow  rate  in  gas  turbine 

c.  generator  field  current 

d.  propeller  speed 

e.  field  motor  current. 

In  the  computer  memory  the  following  constants  will  be 
stored: 

a.  generator  constant  K^ 

b.  motor  constant  K 

m 

c.  constant  of  voltage  difference  N  =  33 

The  required  outputs  depend  on  the  control  system  to  be  used. 
For  an  effective  control  of  the  generator  field  current, 
it  was  decided  that  a  positioning  servo  would  be  used  to  vary 
the  field  generator  resistance. 

D.    DESIGN  CONSIDERATIONS  OF  THE  POSITIONING  CONTROL 
SERVOMECHANISM 

A  positioning  servomechanism  can  be  designed  as  a 

second  order  system  [13  and  14] ,  as  shown  in  Figure  8  with 

an  open  loop  transfer  function 

G(s)  =  R(S)   -    Wn' 


R(s)      s(s+2jwn) 

The  requirements  for  such  a  system  are: 

a.  fast  settling  time 

b.  restricted  overshoot 
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c.   almost  critical  damping. 
Selecting  for  our  system 

settling  time  t   =0.01  sec 

3       s 

maximum  overshoot  M  =10% 

P 

and  J  =  0.9  we  get  from  the  relation 


s   j  .w. 


n 


solving  for 


w     =  444.5       rad/sec 
n 


Therefore   the   open   loop   transfer   function  becomes 


G(s)   = 


R(s) 


197527 


"ClsT         s(s+800) 


The  closed  loop  system  is  shown  in  Figure  8 


HC&) 


<J> 


H759? 
SCs-t  8oo) 


— >■ 


Figure  8.   Closed  loop  positioning  servo 
The  transfer  function  of  the  closed  loop  system  will  be 


R(s)   _      W„ 

■crs5 — rr~r  2 

s  +2iw  +w 
J  n     n 


19.76  10 
s2+800s+19.76   104 


E.         FINAL    REALIZATION   OF   THE   ^-CONTROLLER       ? 

The  microprocessor  based   controller   is   shown   in  Figure 
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Figure  9.   Condition  monitoring  and  control  of  propulsion  plant 

The  output  of  the  system  will  be 

a.  output  1  -  command  to  reduce  fuel  flow  rate 

in  gas  turbine 

b.  output  2  -  value  of  calculated  field  generator 

resistance  for  maximum  deceleration. 
Outputs  are  also  provided  for  alarm  purposes  or  self-test 
diagnostics. 

The  operation  of  the  system  will  be  as  follows: 

Upon  issuing  the  command  CRASH-STOP  the  controller 
is  enabled  for  CRASH-STOP  maneuvers.   The  necessary  field 
generator  current  is  calculated  from  Eq.  (39)  and  then  the 
required  resistance  value  from  the  relation 

Rf  =  3£ 

Rf  -  If 

Simultaneously,  the  W  =  idle  is  issued  to  a  solenoid  valve 
to  cut  the  fuel  flow  rate  in  the  gas  turbine.   Feedback  is 
provided  so  that  a  false  or  slow  response  due  to  false  operation 
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can  be  identified  and  corrected  if  possible.   When  the  pro- 
peller speed  reaches  13  rmp  the  system  returns  to  normal  astern 
maneuvers  by  giving  the  command  to  increase  the  fuel  flow  rate 
in  the  gas  turbine. 

The  above  microprocessor  based  controller  has  been  imple- 
mented for  the  CRASH-STOP  maneuvers  only  but  there  is  no  reason 
why  such  a  system  could  not  be  implemented  for  the  complete 
control  of  the  propulsion  system. 

A  complete  diagram  of  the  closed  loop  system  is  shown 
in  Figure  10.   The  results  of  the  closed  loop  simulation  are 
shown  in  Figures  32  to  35. 

In  particular,  Figure  32  shows  the  variation  of  ship's 
speed  versus  time.   Comparing  this  figure  with  Figure  31 
where  the  open  loop  system  velocity  variation  is  plotted,  it 
can  be  seen  that  the  time  to  achieve  zero  velocity  with  the 
closed  loop  system  has  been  increased  to  24.2  seconds,  i.e., 
an  increase  of  4%  more  time  than  in  the  open  loop  system. 
This  is  to  be  expected  since  the  response  of  the  system  takes 
into  account  the  time  response  of  the  positioning  servo- 
mechanism  involved  for  the  regulation  of  the  generator  field 
current. 

Figure  33  shows  the  variation  of  propeller  speed  versus 
time.   A  comparison  with  the  results  shown  in  Figure  28  for 
the  open  loop  system  shows  a  slight  difference  in  the  time 
required  to  reduce  propeller  speed  from  maximum  to  zero.   To 
be  specific,  approximately  a  2%  increase  in  time  is  required 
for  the  closed  loop  system. 
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Figure  34  shows  the  response  of  the  positioning  servo- 
mechanism  where  R   is  the  commanded  position  from  the  micro- 

c 

controller  and  Rp  is  the  actual  value  of  the  field  circuit 
resistance. 

Finally  Figure  35  gives  the  closed  loop  gas  turbine  speed 
response  where  it  can  be  seen  that  the  increase  in  speed  is 
up  to  3650  rpm,  i.e.,  an  increase  of  1.5%  above  the  rated 
speed  in  contrast  with  2%  increase  in  the  open  loop  system. 
This  is  due  to  the  fact  that  the  small  time  delay  of  the 
positioning  system  causes  a  slight  decrease  in  the  armature 
current. 

An  attempt  to  use  an  operator  controlled  open  loop  system 
for  the  CRASH-STOP  maneuvers  is  considered  highly  inefficient 
since  the  operator's  response  to  fast  changes  in  system  dy- 
namics will  be  very  slow.   On  the  other  hand,  for  normal 
operating  conditions  and  changes  of  propulsion  plant  dynamics 
an  operator  could  quite  effectively  control  the  system  when 
time  response  is  not  a  serious  factor. 

Thus  the  proposed  closed  loop  system  using  a  microproces- 
sor based  controller  is  considered  to  be  the  most  efficient 
way  for  controlling  the  propulsion  plant  dynamics  under  any 
operating  condition,  since  this  controller  could,  as  stated 
earlier,  be  implemented  for  control  of  the  system  under  any 
operating  condition  without  much  greater  effort  for  program- 
ming it. 
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X.   THE  ROLE  OF  SUPERCONDUCTORS 

The  relative  flexibility  and  simplicity  accruing  to  elec- 
trical power  systems,  in  comparison  to  mechanical  systems,  has 
long  been  recognized.   Specifically,  electrical  power   and 
transmission  systems  provide  flexibility  of  installation,  are 
easily  suited  to  automation  and  have  a  high  degree  of  dependa- 
bility [15]. 

Despite  these  advantages,  and  others  that  could  be  obtained, 
the  full  benefit  of  electrical  power  systems  has  not  been 
realized  in  marine  propulsion  applications.   To  be  sure  elec- 
tric drive  systems  have  been  used  in  ship  propulsion,  however, 
such  use  has  been  mostly  limited  to  small  or  intermediate  size 
propulsion  plants.   The  primary  reason  that  electric  drive 
systems  have  not  been  adopted  for  large  power  plants  has  been 
that  the  weight  and  space  for  such  systems  is  greater  than 
that  required  for  geared  propulsion  systems. 

As  mentioned  above,  an  electric  propulsion  system  affords 
greater  flexibility  of  installation.   In  contrast  to  a  conven- 
tional geared  propulsion  system  extensive  shafting  is  not 
required  and  the  location  of  the  power  source  is  not  as  re- 
stricted.  In  addition,  the  large  noisy  reduction  gear  could 
be  eliminated.   The  noise  generated  by  this  component  is  of 
increasing  military  significance. 

For  some  time  now  it  has  been  realized  that  when  the 
temperature  of  an  electrical  conductor  is  reduced  to  a  few 
degrees  absolute,  the  resistance  of  the  conductor  decreases 
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to  an  immeasurable  amount.   The  implications  of  this  fact  have 
considerable  value  for  use  in  electrical  machines. 

If  the  resistance  in  the  windings  of  an  electrical  gener- 
ator or  motor  is  removed,  much  greater  current  densities  can 
be  achieved  with  an  accompanying  increase  in  the  magnetic  flux 
densities  within  the  machine.   If  this  increase  in  flux  density 
is  obtained,  greater  power  can  be  generated  within  a  given 
volume.   Thus,  it  follows  that  the  volume  of  an  electric 
machine  can  be  reduced  if  superconducting  windings  are  used. 

Superconducting  machines  also  have  disadvantages,  how- 
ever.  In  addition  to  the  requirement  that  a  superconductor 
be  maintained  at  a  low  temperature,  there  is  a  limitation  on 
the  amount  of  current  carried  that  is  exposed  to  a  magnetic 
field.   If  the  magnetic  flux  density  increases  then  the  per- 
missible current  density  is  reduced. 

Similarly,  the  very  fact  that  a  superconducting  machine 
develops  extremely  intense  magnetic  fields  invites  yet  another 
complication.   Since  fields  of  the  order  of  several  kilogauss 
can  be  expected  within  the  machine,  disturbance  of  equipment 
in  the  vicinity  of  the  machine  will  certainly  result,  as  well 
as  unbalanced  loads  on  the  superconducting  field  winding. 
Therefore  it  is  necessary  that  a  shield  be  provided  to  confine 
the  magnetic  fields  within  the  machine. 

From  Ref.  15  it  was  found  that  if  a  superconducting  motor 
were  to  be  used  for  the  model  ship  derived  in  this  study, 
certain  variations  on  the  values  already  derived  in  Appendix 
A  would  be  necessary. 
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Table  III  gives  a  comparison  of  the  values  used  on  the 
model  ship  and  those  that  could  be  used  if  superconducting 
motors  were  installed.   The  values  for  the  superconducting 
machinery  were  extrapolated  from  Ref.  15,  Table  3-1.   These 
values  are  the  results  of  a  computer  optimization  method. 
For  a  generator  set  no  data  were  available  from  Ref.  15  but 
it  was  mentioned  that  the  diameter  of  the  machine  decreases 
rapidly  as  speed  increases.   For  the  model  ship  the  generator 
rotor  was  selected  to  be  4.9  feet  in  diameter  which  is  reason- 
able and  compatible  with  that  of  a  superconducting  machine. 

This  study  on  superconducting  machines  was  made  in  order 
to  be  able  to  compare  the  values  of  the  inertia  used  in  the 
model  study  with  another  model. 

It  is  the  author's  belief  that  if  superconducting  machines 
were  to  be  used  instead  of  the  conventional  DC  machines,  then 
the  variation  of  the  specified  parameters  in  Appendix  A  would 
be  minimal  and  the  results  would  not  differ  much  from  those 
obtained  in  this  thesis. 

If,  on  the  other  hand,  conventional  DC  machines  were  to 
be  used,  modification  of  the  calculated  parameters  most  probably 
would  be  necessary,  since  all  calculations  were  performed  on 
data  extrapolated  from  different  machines. 

It  appears  that  the  results  obtained  in  this  thesis  for 
a  conventional  electric  propulsion  system,  would  remain  sub- 
stantially correct  if  superconducting  machines  were  used 
instead. 
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Table  III.   Comparison  Between  Conventional  and  Superconducting 
Machines 


Operating  voltage 
Operating  current 
Motor  speed 
Approximate  diameter 
Weight 


SUPERCONDUCTING 


300  volts 
50,000  amps 
200  rpm 
6.39  ft 
113,000  lbs 


CONVENTIONAL 
MACHINE  USED  IN 
MODEL 


400  volts 
33,000  amps 
2  30  rpm 
6.5  ft 
90,000  lbs 
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XI.   CONCLUSION 

This  thesis  has  presented  a  control  model  for  dynamic 
braking  of  ships. 

The  use  of  different  braking  devices,  such  as  a  resistor 
bank  or  a  mechanical  brake  has  been  avoided  by  the  use  of  the 
inertia  of  moving  parts  of  the  gas-turbine-generator  set  as 
load  on  the  motor,  which  during  the  dynamic  braking  phase 
effectively  operates  as  a  generator. 

It  has  also  been  shown  that  by  effective  control  of  the 
generator  field  current,  maximum  deceleration  of  the  propeller 
shaft  can  be  achieved  without  exceeding  the  overspeeding  limits 
of  the  gas  turbine. 

The  model  used  permits  the  study  of  different  motors  and 
generators  as  long  as  their  characteristics  are  specified. 
It  also  permits  the  evaluation  of  a  given  propulsion  plant 
with  different  types  of  propellers  if  their  characteristics 
are  also  known. 
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XII.   RECOMMENDATIONS 

The  control  method  that  has  been  presented  for  the  auto- 
matic control  of  ship  deceleration  during  emergency  conditions 
provides  a  realistic  tool  for  further  studies.   The  following 
are  topics  suitable  for  future  investigation: 

1.  Design  of  a  microprocessor  based  controller  for  all 
phases  of  propulsion. 

2.  Effect  of  sea-state  conditions  on  the  operation  of 
the  propulsion  plant  during  CRASH-ASTERN  maneuvers. 

3.  Development  of  a  complete  propulsion  system,  in- 
cluding gas  turbine  dynamics  and  control  based  on  microprocessor 
based  controllers. 

4.  Utilization  of  superconducting  electric  machines  in 
propulsion  systems. 
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Fig.  11-Valucs  of  the  theoretical  braking  resistor 
for  maximum  deceleration. 
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Fig.l2-Approximated  values  of  the  braking  resistor 
for  maximum  deceleration. 
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Fig.lU-  Plot ■ of  the  propeller  speed,  vs. time  during  dynamic 
braking  using  the  theoretical  and  the  linearized 
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Fig.l5-Power  dissipation  diagram-Theoretical  resistor 


PMAX=Maximum  power  dissipated 

PRB  =Power  dissipated  in  braking  resistor  RB 
PRM  =Power  dissipated  in  internal  resistor  RM 
PERCl=Percentage  of  dissipated  power  in  RB 
PERC2=percentage  of  dissipated  power  in  RM 
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IB-Power  dissipation  diagram-Linear  resistor 


PMAX=Maximum  pcser  dissipated 

PRB  =Power  dissipated  in  braking  resistor  RB 
PRM  =Power  dissipated  in  internal  resistance  RM 
PERCl=Percentage  of  dissipated  power  in  RB 
PERC2=percentage  of  dissipated  power  in  RM 
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FiR  17-Plot  of  ship"s  velocity  when  a  fixed  resistor 
g     is  used  for  dynamic  braking. 
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Fig.l9-Power  dissipation  diagram-Fixed  resistor. 


PMAX=Maximum  dissipated  power 

PRB  =Power  dissipated  in  braking  resistor  RB 
PRM  =Power  dissipated  in  internal  resistor  RM 
PERCl=Percentage  of  dissipated  power  in  RB 
PERC2=Percentage  of  dissipated  power  in  RM 
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Fig.21-Plot   of   propeller    speed   vs. time   when   an   air 
brake    is    used    for   dynamic    braking. 
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Fig  22-Plot  of  propeller  speed  vs. time  when  a  combi- 
nation of  a  variable  resistor  and  an  air  brake 
is  used  for  dynamic  braking. 
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FiS-   26  -Plot  of  ships  velocity  vs. time  when  a  combination 
of  an  air  brake  and  a  fixed  resistor  is  used  for 
dynamic  braking 
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Fl^.  25.      -Plot  of   propeller   rpm   s  vs. time  when  the   inertia  of  moving 
parts   of   the   gas-turbine-generator   set  is  used     as   load 
for  dynamic   braking 
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Fift.  30-    -   Plot  of   gas-turbine   rpm   s  vs. time  when   the   inertia   of  moving 
parts   of   the   gas   turbine   generator   set   is  used   for   load 
during  dynamic   braking 
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Fi$.  31-    (-Plot  of   ships  velocity  vs. time  when   the   inertia  of  moving 
I  parts  of   the   gas   turbine-generator  set   is  used  as   lodd 
during  dynamic   braking 
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F»^.   32.     -Closed  Loop-Plot  of   ships  velocity  vs. time  during  dyruhnic 
braking  and   astern  motion 
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Fig. 34-Closed    loop   response-Variation   of    generator 
field   resistance 
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Fig.36-Closed    loop   response-Variation   of   the    generator 
field    current    vs. time. 
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APPENDIX  A 

DETERMINATION  OF  ELECTRICAL  AND  MECHANICAL 
CHARACTERISTICS  OF  THE  SYSTEM 

A.    DETERMINATION  OF  THE  EL6TRICAL  CHARACTERISTICS  OF 
ELECTRIC  MOTOR  AND  GENERATOR 

Since  data  of  an  actual  working  system  were  not  available 
at  the  time  this  study  was  conducted,  an  effort  is  made  here 
to  establish  the  electrical  characteristics  of  the  motor  and 
generator,  based  on  general  known  characteristics  [4]. 

The  maximum  torque  on  the  propeller  due  to  the  water 
moving  through  its  blades  can  be  found  from  the  tabulated 
data  of  the  propeller  and  the  equation 

Q  =  C   .p.  D3  (V  2+  (N.D)2  )  (1? 

vp   q  v  p 

Thus    for    the   maximum  number   of    shaft   rpm's    this   was   "found    to 
be   equal   to 

Q     =583580        ft-lb  (2) 

P 

Assuming  the  frictional  torque  on  the  shaft  to  be  constant 
with  a  mean  value  of  22,500  ft  lbs,  then  the  maximum  torque 
exerted  by  the  motor  according  to  the  equation 

2TTI  Up  =Qm  -  Qp  -Qfr  (3) 

under   steady  conditions    is    found   to  be 

Qm  =  Qp  +  Qfr     =   608580     ft-lb  (4) 


99 


Translating  this  to  horsepower  gives 


HP  =    ■  g3jj  ?     =   26650.87     Hp  (5) 


Expressing   this   value   in  kilowatts  we   get 

P  =  Hp  x   .746   =  19881.55     Kw  (6) 

For  marine  applications  the  maximum  operating  voltage  of 
motors  is  600  V.   Therefore  the  circulating  maximum  armature 
current  will  be 


P( watts) 


=  -777 — -,,    <t      =   33135.9     amperes  (7) 


By   now  motor   efficiency   of  1    has   been   assumed. 

For   this    size   motor   almost   6%   of   the   power   is   dissipated 

in  the   armature   resistance,    so 

Rm  +  Rg  =  R     =      33qqq00      =  0.00109     ohms  (8) 

Assuming  equal    internal   resistance    for   generator   and  motor 

we  get 

Rm  =  Rg  =   5  10"4  ohms  (9) 

The   back   emf  of   the   generator   becomes 

Eg  =   600  +  la  .   Rg  =   616.5       Volts  (10) 

and   of   the   motor 

Em  =   600  -  la   .   Rm  =   583.5         Volts  (10a) 

under   steady   state   conditions. 

So   for   ahead  maneuvers   at  maximum  speed 
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Eg  =  Em  +   33       Volts  (11) 

and   for   astern  maneuvers 

Eg  =  Em  -  33         Volts  (Ha) 

The  generator  back  emf  on  a  D.C.  machine  is  given  by 

E  =  k  .if  .  n  (12) 

from  which  solving  for  K 

K  =  E  /(if.n)  (13) 

Assuming,  at  steady  state,  a  field  current  of  10%  -  12%  of 

the  armature  current 

then     Ijr     =  40   amperes 
r  max 

Therefore  at  maximum  operating  conditions  at  steady  state 

Em  =  583  volts 

N  J  =  60  rp/sec 

gt 

N     =  230  rpm 
prop        c 

Eg  =  616  volts 

So 

Km  =  0.06326    and    Kg  =  0.25666 

Similarly,  from  the  torque  equation  of  a  D.C.  machine 

Q  =  K*  .if  .ia  (14) 

having  in  mind  that  the  output  torque  of  the  motor  must  equal 
the  sum  of  friction  and  propeller  torques,  is  found  to  be 

Q_       =  608580     ft-lb 
TTiax 
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Solving  Eq.     (14)    gives 


Km     =  0.459153 


For  the  generator,  considering  maximum  operating  conditions 
and  steady  state  at  3600  rpm  the  tabulated  data  give 

Qg  =  38000   ft-lb 

therefore    from  Eq.     (14) 

•Kg'    =  0.0286697 

The  above  results  are  summarized  in  Table  IV. 

B.    CALCULATION  OF  PROPELLER  MOMENT  OF  INERTIA 
From  [16],  propeller  weight  is  given  by 

W  =K.D3(MWR)(BTF) 

where  D  is  the  diameter  of  propeller  in  inches 

D  =13.94  ft  =167.28  " 

MWR  is  mean  width  ratio  =  .8  for  a  three  bladed  propeller 
BRI  is  blade  thickness  fraction  and  is  assumed  to  be  .05 
K  is  assumed  to  be  .26  for  three  bladed  propellers. 

Substitution  yields 

W  =  48681.66     lb 
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Table  IV.   Summary  of  Characteristics  of  Motor  and  Generator 


AA 


DESCRIPTION 


MOTOR 


GENERATOR 


1  Torque  constant 

2  Electrical  constant 

3  Back  emf  (max) 

4  Internal  resistance 

5  Internal  inductance 

6  Field  voltage 

7  Field  resistance 

8  Field  inductance 


Km  =0.459153 
Km=0. 06326 


KG  =0.0286697 
KG=0. 0042916 


Em=583.5  volts    Eg=616.5  volts 


-4 
Rm=5  10    ohms 


negligible 


E,  =400  volts 
fm 

R  £=10  ohms  min 
mf 

L   =20  hrs 
mf 


-4 
Rg=5  10    ohms 

negligible 

E_  =400  volts 

fg 

R_  =10  ohms  min 

fg 

L,  =20  hrs 

fg 


At  steady   maximum  field  current 
state 

maximum  armature  current 


I«-  =  40  amps 
I   =  3  3136  amps 


Notice  that  the  difference  in  Eg  and  Em  is  3  3  volts;  therefore, 
without  overloading  the  armature  circuit,  application  of  astern 
motion  starts  when  Em  =  33  volts. 
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Therefore  propeller  mass 


m=—=  1511.85   lb-sec2/ft 
g 


The  radius  of  gyration  is  assumed  to  be  .25  D  so 

I   =  (.25D)2  m  =  18361.7  ft-lb  sec2 
P 


An  allowance  of  25%  increase  is  allowed  for  the  inertia  of 

entrained  water  in  the  propeller.   Therefore 

1=1  (1+0.25)  =  22952.2  ft-lb  sec2 
P    P 


C.  COMPUTATION  OF  MOMENT  OF  INERTIA  OF  MOTOR  ROTOR 

Since  no  data  were  available  for  a  D.C.  motor  this  size, 

considering  data  which  apply  for  a  synchronous  motor  of  15000 

Hp  rating   output  power  and  which  give  the  weight  of  the 

rotor  to  be  approximately  69,000  lbs,  a  rotor  weight  of  90,000 

lbs  was  selected  [12] .  Therefore 

m  =  —  =2795.03  lb  sec2/ft 
g 

and 

I  =1/2  m  r  2 
m        I 

For  an  outside  diameter  of  the  rotor  of  (  2  m  )  6 . 5  ft  we  get 

I  =  15042.72   ft-lb  sec2 
m 

D.  COMPUTATION  OF  MOMENT  OF  INERTIA  OF  PROPELLER  SHAFT 
Assuming  a  shaft  length  of  50  ft  with  outside  diameter 

D2  =  1.6  ft  and  inside  diameter  Dl  =  1.37  ft  gives:  [12] 

Volume  =  2n(r2-R1)L  =36.12   ft3 
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The  specific  weight  of  steel  is  490  lb/ft 

so  weight  of   a   shaft  w  =  17702       ^ 

m  =   549     lb 
so  mass 

and  I  =   175.9     ft-lb  sec2 

Total   moment   of    inertia   of  motor-propeller    system 

I  =  Ip  +  Im  +  Is  =   38171         ft-lb  sec2 


E.  CALCULATION    OF    GENERATOR    ROTOR   MOMENT    OF    INERTIA 

As  in  the  case  of  the  motor  and  taking  into  account  the 

maximum  number  of  revolutions  of  the  generator  a  rotor  weight 

of  40,000  lbs  [12]  was  assumed 

Therefore 

m  =  —  =1242.24   lb 
g 

and  assuming  an  outside  diameter  of  D  =  4.9  ft  we  get 

Ig  =  1/2  m  R2  =  3760.7    ft-lb  sec2 

F.  CALCULATION  OF  TURBINE  ROTOR  INERTIA 
From  [9]  the  following  data  were  extracted 

Low  pressure  compressor  Jp  =  586 
High  pressure  compressor  Jp  =  4  89 
Free  turbine  rotor  Jp  =  5009 

where  Jp  =  polar  moment  of  inertia.  Therefore 

T  ID  5009  1CC     rc  ru.    1K    _       2 

1     =12.2  =^0  =  155'56    ■   ft-^sec 
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Therefore  total  moment  of  inertia  of  gas  turbine-generator 

system 

Igt  =  Ig  +  Igen  =  3916.243    ft-lb  sec2 

G.    CALCULATION  OF  THE  WEIGHT  OF  BRAKING  RESISTORS 

In  order  to  be  able  to  select  the  system  that  will  be 

applied  some  criteria  about  the  size  of  the  braking  resistors 

are  needed  for  each  case  [17] . 

The  energy  that  has  to  be  dissipated  in  the  resistor  is 

given  by 

2 
E  =  i  (t)  xR,     in  joules/sec 

The  energy  .which  has  to  be  absorbed  is  given  by  the  area  under 
the  curve  of  the  previous  equation. 

If  the  resistor  were  of  copper  with  specific  heat  of 
0.0918  BTU/lb°F,  for  an  acceptable  temperature  rise  of  100 °C 
the  weight  of  the  resistor  will  be  given  by  the  formula 


„  .  ,  .  .  ,,     Energy  absorbed  in  joules 
Weight  in  lbs  =  "  RTTI J 


0.0918- 


BTU 


^x212  'F  x  1054.8  l~?s 


ih,-r      *««-r.«  BTU 

The  energy  absorbed  in  joules  by  the  braking  resistor  in 
each  case  was  found  by  using  Simson's  discrete  integration 
formula . 
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Table  V.   Weights  of  Necessary  Braking  Resistors 

Method  of  Energy  Absorbed     Weight  of 

Dynamic  Braking  in  Joules  Resistors  in  lbs 

Use  of  a  Variable        1.38717  10  675.7 

Resistor 

Use  of  a  Fixed  1.7621  106  85.883 

Resistor 

Use  of  a  Mechanical       1.3870  10  673.3 

Brake  with  a  Variable 

Resistor 

Use  of  a  Mechanical      1.2549  106  61.131 

Brake  with  a  Fixed 

Resistor 
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#. -Map  of  propeller  torque  characteristics 
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Fig. 39-Propeller  torque  characteristics ,Cq 
as  a  function  of  the  second  modified 
advance  coef f icient  ,s 
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H.-  Map  of  propeller  thrust  characteristics 


Fig. 40-Locus  of  thrust  coefficient  Ct  and  second  mo- 


dified advance  coefficient. 
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5.-  Map  of  ship's  resistance , propulsive 

coefficient  anc  propeller  speed  vs. 
shJ  p  '  s  spe*-C  . 
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Fig  m-Plc".  :■-'  S-.x^'s  resistance  (R ),  propeller 
speed  (Np)  and  propulsive  coefficient 
(n  )  versus  ship's  speed  (V). 


112 


APPENDIX  B  -  COMPUTER  PROGRAMS  USED 


v  ThIS  PROGRAM  »vAS  WRITTEN  IN  nsi/360  SIMULATION 

*         LANGUAGE  .  TRANSLATION  Of-  THE  PROGRAM  WOUI  G  ^e 

1 Pr')'JTPj:c   Tr   anothtp — r,i."ut  at  icn    l  AMr.iiar.F    Ts    i  sfp 


TITLE.  FRCGRAM  A 


*         rpFM  LGCP  SYSTEM 

i       

-* nr-Tfc-piU'NATlur"    CT-     THE    VALUES    Q-f     THF     TWrij^PTTCAt 


*  VARIABLE    RESISTOR 

■v 


V 


*  VARIOUS    CASES    OF    DYNAMIC    BRAKING    ARE    INCLUDED    IN    Th'^ 

*  PROGRAM     ,     AND    ARE    THE    FOLLOWING     : 


1. -DYNAMIC     BRAKING    USING     A    LINEAR    VAP*AaLF 


* 

t rf 3*  srrrft 

*  2. -DYNAMIC  BRAKING  USING  A  FIXED  RESISTOR 

* 3.-rmiA-',T,    tpaktng   u;;i\r   a   "echan'Cal 

~  a;„     tYPl-     bPAK.E     . 

*  4. -DYNAMIC    BkAKiNG    USING    A    COMBINATION 

-a U£ — A    VAST  ARE  E    '■'  I  S  i.STQR — AND     THC     MEGHAINQ  (",  M 

*  -  ATP     BRAKE 

*  5.-DYNAMTC    BRAMNT.    USING     A.    COMBINATION 

p rf-e — ifti — (Vfft — B  PAKE — ANP    A — F  '  X  E  o    P  F  S I  S  T  P  R — ; 

"*  '.,HAT     TEE    USER     HAS    Tfj    fTD     j  S    TO     SELECT     A     METHPO     RY 


Frp     EXAMPLE- G  A  S  j  =  2 M  i  A  N  S    T  M  A  T    0  y  N  A  ^  T  C — bRAK,  T  N,  r 


WTi  L    HE     AC.rO^PLTShEn    BY     THE    USE    OF    A    FIXED 
P  E  ST  STOP 


IN    ADDITION 


*  THE    USE'-    HAS     TO    SPECTnY    The    FCLLOWING     : 

■ (_A_) m yAX'r'll"     PFWOTSSIRLfc     CURRENT .N     TRc-     ARMATURE 

*  C':RCLtT   (THIS  :s  only  useh  in  the  calcu- 

*  LATION  HP  THE  THECRCT\CAL  VALUES  OF  THE  VAST  ABLE 

I krsrsTCR    ) 

*  (B)     -    VALUES    Or    THE    ELECTRICAL    CONSTANTS     CF     THE 

*  vpTC     AND    GF^FRATOR.hM     Ai\C    KG     c  E  S  PE  CT  j  V  EL  Y  .  

I rn    -    VALU"^    OT     THE    FTEXTTATTOi:    CONS'  PJT3    UP    M'jTijP 

*  AND    GENERATOR. KMM    ANT    KGG    R ESPECT [ Vp L Y     . 

■ UCU -     VALUED     OF     THE — INTERNAL    k  ES  "  S  TANG  E  S     CF     "rITPK 

*  ANn    GENERATOR*    RM    ANC    RG 

*  (E)    -    VALUES    FOR    THE    INERTIA    C^    PROPELLER     , 

B SHA^T     frfte    MOTOR     SYSTEM     AND VALUES     ri> — Wf? 

*  INERTIA    OT    GENERATOR     RCTCR     ,POWEP    TURBINE     SVSTEM 

*  AS    "N     APPFNO.-X     A    CF    THE    THESIS; 


V 


-iNTLG    TRAP? — 
iiSTGER    NFLtJT 
tNTGE R     CAS F 


"LtN^r    i\Plf)T=l 

-i ^f-TcpMfVATTy-I.M     Q£     P»np-|QRE     PARAMETERS 


CLNST     v=?-£il  71 .  ,C  =  13.<34 
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,CNS1 


nFJ^RMJNATION    OF    SH7P    PAP  a M E T E P S 
! .  7  3  9  7  E  +  C  5 


rETERMUATION    OF    MOTOR    CHARACTERISTICS 

,Li\  -jT    K»i  -  0  *  H.3  P.C  t  v.  >***  0  ,A  b  0 Lr:  2  ,  R M  -  6  .  GO  0  5 


DEJTRJMT_MJ_iofj    pc       GENERATOR    Ch ARACTERISTICS 


,CNST    KC  =  rJ.?  5  6  6  6^  ,  KGG=0  .028669  7  ,  RG  =  C  .00C5 


iNFRTTA    GF    GAS-TURBINE/ GENERATCR    SET 
rfiST     .■.C-1nl6. . 


CTHFR    CFNSTANTS 


GNST    P f=2.14l5  9,P=?. 


i N 7 T * a U    C  C  m~  !  t  : 


!NS 


NLCN     'Ci=30, 


!  C 2=3.8333 1 1 C3=0.0, 1 C 4=60. 


ER'IVAT.TVE 

l  =  'm-;qi  LLCJ  lilllCJJ 


K=INTGRL  (TC3,UJ 

fJGT  =  INTGRL  (  IC4,GT0C)T] 

P  T  -  t  p  3  -i  (  v  ) 


b=TA6W(  V  ) 

VP=  (  1. .  -W  )  *v 

r  =  i  Ai-.r  (  v) 

ft=VP**?  +  (N*C  )  *  *  2 
aaS-QBJLLAJ t 


CT=TABCT(S,VP) 

T  a  -  r,  t  ■-  p  a  ( i-  a  a  g  ;  a  a 


TP=(1.-T)*TA 

tpr=itp-rt)/m 

utct  =  tijr 

CG=TABGQJ  SjjVP) 


R=INTGRL (IC2.ND0T) 

|  ,v  -  :<  (v  a  T  f  v  *  f\  p,  f  a  i 


EG=KG*IFG1*NGT 

(Ar"325c 

rf  c  1-  [r"-:-;. )/(  Kr&*w?f ) 

q»c:=-~  a*  7  c  r/^KWM 

P^PMIA**2J*RB 

P  p  v  3  ( T  A  5  *  ?  1 I  *  ft  jvj 


PVAX=CM*IA 

yf|_T-!-M-T  AtOV 


PFRCl=( PRB/FVAX )  *1  00. 
PFRCZ^ ( PRM/FMAX»*100. 

•  T  *  >'»'  -  T  '.'  W  g : 


n\AiV{c 


■('J.LT.0.216u)     &6    TH    700 


TF(rASF.Fr;.l  )    GO    TO     100 
IHCASF.Eq.2)    GJ    TO    200 


MUStr.bC.'i)    GU    TO    TUTJ 
IF(CASE.EQ.A)    GO    TD    200 

-f  F(  TASF  .  FC  .5  )    GjO     TH     330 
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100 


j FV=40. 


,E*=( EM-IA*RN )/! A 
Flu;  -i  T.r.nms  )    Sfl    to    ]  in 


Gr  TO  1.2C 
11C  RB=3.0315 
t£0   cc:c:t  rMUF 


SG    Tn    600 


|Sr''MMPr.T     '"ATM     nrNFRATPk     AND     ALR     A    CONSTANT     FLS'STQP 


iCO    RP=0.017 

i  r  c  =  ■>  o » 


i.vcT=o.n 

SC    TO    6  0C 


St   CF 
J00    P< 


AN 


9  1  » 

5.* 
=45 


TP     CLUTCH    eBAKING    DEVICE 


CL 


TIM 

50.* ( FS-PC) 


im:t  =  3.o 


:f(ccl.lt.c.o   go  to  301 
;r   in   30? 


201 
202 


CL 


-0. 


0 
U6 


P  IV  t-I  T  - 1  5  0  ■ 


r 

2G3    (, 
iM   C, 


CL 
F( 

t~ 
CL 

flKI 


CCl 


55 
G 


Trr- 
=  QC 
T  TM 


b.=>FNFT 

T.CCLM)     GO    TO    303 

n 


TV. 
i  M 

Of 


F(C/^E.FC4)    GO    TC    100 
;r-(CASE.EC5)    GC    TO    200     • 
iC    V)    600 


fcOC    Sl'H1  =  'JMH,-CP-OI-P.-'j(,l+gMt+uM(;b 

GO    TO    900 
70C    7F<C^SF.GE.3.  AMC.CAS1".  1,  T  .  fr  )     CC    TO 


OU- 


QCL=0.0 
CC    TP     711 


10    PC-r.  K).,?124*(N**Z) 
K=TTME-6.7 
PSI  =  P  S - 3 G . * K 


U',  ^SMt.'M^I-Pf.  J 

IF(CCL.LT.C.O)     GO    TO 

1£L    7j 


'0  1 


GC    Tf 


^L2- 


7C1    GCL  =  n.o 
702    CrNT  i!MUE 


-ill  C'Vr,-0.  T 
PP3=0.0 
FRV=0.0 


"PTPTT 
PFPC2- 


0.0 


EM=0.0 
RB=C0 


i/>ST=-33O0C. 
toA^T-KMf>*-F^*IAS  [ 


8C0    SIJMO=QAST-CP-CFP- 

Uip  ccr^-^niF 


•cci 
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CPF  =  f-CNS!rf(N  .-22500. 

|\ r n T=C>U^C/  P.^I^  ) 


0.0, +22500. 1 


T^ 


r*r- 


j)pfn=poB/  iconoo. 


p^A  <L  =  ^AX/10JOO  ), 
li\\=:  A*C.C1  515 
iur,\m  j-ir,/\00OO, 


cp2=cp/ic-:cc. 
el  1=0cl/icc0c0. 


Li_ 


KL    F!NnM  =  25.,CELT=.l 


T   .1 


PRBtVCH.TtGPtVf 


net  ■=;  = .  L 

!  A,PRN  ,PFRC1  ,CFR,NF 


AL,RB,PMAX, PERC2»CMG, 

t  ,'.  AST,  OAST  tQMGOtSU''>3,S,  ECCJ,  PS  ,«CL  ,M3  ,J1  F,CT,PC.Qr,PN  ,\GT  ,  . 
lFni,NF,IA,QE»VP»W,T,A,B,TPR,K,PSI 


THE    FOLLOWING    STATE^FMTS    ARE     APPLTf.A:nF    TH     THF 
PI  1TT'^;r:    PACKAGE    FHP    OSL/360     S  .' *'JLa  T  ■  CO    l_  ANgjJ  AG  E 


AT  N.P.S 
OTHER  THAN 


The  STATEMENTS  SHOULD  BE  CONVE! 
THIS  PLOTTING  PACKAGES  ARE  AVA' 


TFD 
LABI 


Cfli.u    riHWr.  (  I  ,1.  ,Tr-*F,V) 
CAlL    OP,WO( 7 ,1 ,TI ^E,MRFAl   ) 


T — r» 


-rr 


RTfffl  *  . 1  pn  TETPTrg 
L     )RWG(  7  ,?,TI,,J,F,  PRM1  ) 


CAI. 


-H- 


CAL 

CAL 

-tAi 


')  R  W  G  {  1  ♦  4  ,  T  I  P  F  ,  P  c  P  C  7 
DPWG{?,5,TIMEt  PMAX1  ) 

J — DPWG(  ■'!  ,  I,  TIME,  i  AT 


-J- 


C£L 

CAl 

CAL 


L    ORV*r<  V,  ?,T  !  ME,^) 
L    DRWG( 5,1, TIME, QCL11 


ENORVI  (\PLOT) 


ttCTt 


TABULATION    OF    SPEED    REDUCTION    COEFFICIENT    VS.CDFED 


FUN 


CT" 


ON    T 

A  y/T/o. 

A    ->IT/4* 


A  H  b  (  V  ) 

\iM13)  ,'.'JT(  13) 


GAT 
OAT 


,5., 7. 5, 10. ,12. 5, 1 5., 17. 5, 20., 22. 5 ,25. ,2  7.5,20./ 


.JG3,.J1,.><Z,.C3  3,.C4C,.0  45,.0:S,.02,.0'~4/ 


JFTTTTJT . 

£  p  |    \j  _  T       e; 

N  =  'FTX(  V/'T^LVi+l 


.  )    GO     Tf|    I 

EJ 


SI.  G 

TAB 


PW=(WT( 
fo=SLCPW 


CT!)f-"'i 


N+l)-WT(N)J/(VT(N+l)-VT(N)) 
*  (V-VT (N) ) +WT(N) 


T  A  !J; 

RET 

TNTT 


w=0 . 0 

UP  N 


TABULATION    OF       SFIIP'S    RESISTANCE     VS. SPEED 


_Flll'CTTriN    TjERlVj 

TJiW-f  v\(!\    VTTTTJ  ,k'!T(l3l 
DATA    VT/0.C,2.5,5.,7.5,  10. 
-i.ATA    pTT/2tQ.Q,  7000.  .  L3Q0JJ 


12. 5, 15. ,17. 5, 20. ,22. 5, 25. ,27. 5, 30./ 

r^nnr.,;cniin.t^7f)nn.,fl[-innr.,if-riiP. 


H3ft.10  3.  ,172C03.  ,2250; 
!F(W.GT.3C.)     GC    TO    I 


J. ,23  3000./ 
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n  f  l  / = ? .  s 


rr 


|ABR=S10PR*  (V-VT  IN)  J+RTTIN  ) 

|£I1LBJ 


(  !'i  )  }  /  (  V  T  (  N  + 1  )  -  V  T  (  M  )  ) 


TAt<R  =  2800CC 

Ifturm 

§£&£ 


TABiJL  AT  TCTi    OF     THRUST    REDUCTION    COEFFICIENT    VS.SPEbC 


FU:CTTON    TART(V) 
D:>'rr!<r."OM     VT  (13)  ,  TTt  13) 


CA 

CA 


N  = 

IF 

31 

TA 


,5.,7.5,10.,12.5,15.tl7,5t20.,22.5,25.t27.5t30./ 
C.0,.01t.05,.07».08».075,-072,.06,.02/ 


TA    7T/0, 
TA     TT/5  = 

i.y-g.g 
T  c  y  x  (  v  / 

( v.r>T.^ 

BT=SLGPT* (V-VT (N) )+TT(N) 


CEl.VJ  +  l 

(\+l    -Tt(hh/ WT(i\  +  l)-vT(M)) 


TA 

RF 


BT=0.0 
TURN 
fl 


TTB1TCFTTDT3  of    riiRusi    ciefi-h;  if. nt    vs.5[_CH^r    mi;i  'i-?tr 

ADVAMCF  CnFFF  ITTFN'T 


fUNCTTON  TABCTl  S  ,  VP  ) 

oh^i;s:gn   st  <2i )  .cTirn  ,trz  1 211 

L1ATA    ST/-l,,-.9,-.8f-.7t-.6t-.5f-.4, 


A 


.l»O.0t.lt. 2*. 3 

It  .4i 


\Frfa  CT^;-.4,-.io,-.  35,.  T5  t.l,.l5,.?,.26,  .36,  .44,  .30,  .: 
1. 45,. 42,. 3 9 t. 37,. 32 1. 33 i. 50/ 

LATA  Cri/-.^t-.2,-..?H,-.;?,-.j5,-.3a,-.4,-.41,-.j6,-.?' 
i- . 3  5 ,  -  .  3 4  ,  -  .2 5  ,  -  .  2  , -  .  1 3  ,  -  .  G5  »  .  0 1 1  .  iO  ♦  .  2  1 »  .  45/ 

CELS=0.i 

i\-r  r?\  (  (  5t  i.  )/ncLr>)  rl 

IF(VP.LT.O.O)    GO    TO    2 

SLrCTl=(CTl(N-H)-CTKN))/(ST(N-H)-ST(M)) 

|Ah(J  =  b!_LC  I  1*  I  S  —  S  *  IM)  )+(.  '  1  IN) 

GC    TP    6 

Si.  rr.  T  2  =  (  C  T  2  ( N  +  l  )  -C  T  2  t  i\i  )  )  I  (  S  T  ( fs  -H  )  -  S  T  (  M  )  ) 


-« «-V *- 


|"ABCT=SL0CT2*(  S-ST(M)  I+CT2  (Ml 
RETURN 

-f*e — 


lABULAi'OT'    'IF     I'JkUUb    CTJFFFICTETTI     ^rrTTXrT"*T~~^_TTF-'TT~ 
ADVANCE    COEFFICIENT 
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■ND 


TABULATION    DF    PARAMETERS    OF    THF    PHI^    vrvE 


HN    TM<06{  V3.WF  J 
iTru    V?T (5 ) , WFT( 10) ,OET( 100*100) 


CC. 
CC. 


=  4 

—■a 


OIL 


=  5" 


GC 
OC 


-  fe  3  C  C  . 


.o: 


Tr- 
ee 

en, 


)  =  12 
1 

_2_ 


SCO. 
0000. 

£  c,  G  C . 


7C00. 

5  C  C  0 . 


~T2r^2, 


h 
7 


=  6 
=  3 


2C00. 

ccoo. 


r£->- 


7C00. 
£500. 


25CC. 

cooo. 


4, 
"5 


'4500. 
■\  500. 


eoco. 

400C. 


1  =  6 
10  )  = 


4CC0. 
6^CC0 


ICC. 

3CCO. 

i-500. 


■>  • 
3 


4 
5 


40CO. 

<3000. 

tCCQ. 


7 


C5CC. 
5C0O. 


-  3  C  C  C  C 


)- 

=  5 


C50CO 
CCC. 


■+  < 


=  1 


:  c  o . 

•coo. 


2G00. 
6  C  0  0 . 


^coo'. 

H500. 


■■4  3000 
ICCO. 


CCCO. 
"030. 


6 
7 


'CCCO, 

:4CC0. 


IT 

9 


LQ-jaZA 


7TCU 

C5C0. 

■>&.rmc 


;F( 
IF( 


V? 


IT 

gt 


.0)    GO    TO    I 

CCC.)     GC    TO    2 
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loo 


X 

G( 
T: 

GA 
Ji 


TO    100 

TO 

to  mo 

LQUL. 


;rj   to 

[  F  ( in  F 

.re..'- 


iFl 

;c 


i  o  o 
.  t  .  3 : 


oo. J   go  to  : 
3-g  j.  )   on   th 


c 


TO 

TD 

=12 


7  00 


FIX 


?00 
20JL- 


(  y 
(  (w 


■  v;t(  • ) 


7TC00TT+X 
F1-330Q.)  /1 000.  )  +  .1 

x_ 


V  i  = 
FL 


V-  F I 
0R= 


r,n  n- 


-W  F 
<  OP 

-rrro 


T(  J) 
/1000.  ) 


(  Q  F  T  ( 
;r;ET<  1 


TAHOE=OET(  t  ,  J)  i-QELQR+OELQw 

K^T||P!fM 


+1, J  )-QET< 
H-l)-:bT  ( 


I, J)) 

.  ,J1) 


END 


THE     rOLLCWNG    DATA    APE     PART    OF    THE    PLOTTING    PACKAGF 
FOR    DSL/360    SIMULATION    LANGUAGE    AT    N    P    S    .     THEY 
'■Iv'-'ij;  r, — pf    "■*on7r>ro    'l1    HTMEP — Tffrft    th"S    PtrftTTWfi    ^r^np 
ts    AVAILABLE 


fi.c.eytn    or   a 


3T    CF    THE    VELOCITY    VS. TIME 
FLANT'NTS  I A 


C er: Ti T-. 5i 

:T  CF  THC  VAQT^TTGN  CF  PROPELLER  SPEEC  VS. TIME 

i-LM-TU;"  S I A 


5T7  -230.  100. 

T    EE      PR^l.PR/vi  ,PERCi»PERC2»PMAX 

'LAM' 


JUHO. 


-4A- 


j-  0.34  0.0 

T    tiE       :  A     AMD    FN     VS  .TIME 


30. 


LAfUf.'J     r, tA— 

J                          I. 

: r  c f  ha  v s . t I v e 

J.O 

15  J. 

5. 

5. 

PLANT7 Nv S              lA 
j                     )  .  3  8 

0.0 

0.004 

5. 

5  . 

— 
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fcXtC    r)r-L 
n SL.  :NPUT 


in    * 


HIS 

A  NO 


PRCGRAM     I'JAS     WRITTEN     IN    DSL/3o( 
IAGE     .     TRANSLATION    CF     THE     FROG! 


E     ANGTITn — SIMLLATVGN     L  AN OH  AGE 


STMiJl  ATI  ON 
,M     WOULD    t^E 


'S     IJ  SF  0 


p  r  r  r-  p  a  i   b 


f.P^M     I  TP 


THF  FOLLCWING  PROTRAM  CALCULATES  THf1  OPEN  LCOP 
-tW- TuE  SYSTEM  UNDER,  U  Y  "■!  A  ■  .'IT  C  HRAKlfvC  C  C  Ml  I  T  I  ON  S  . 


RFSPCNS.f 


A  L 
CRG 


I  N  E  A  R 
M  THF 


RR.AKT!\G  RESISTOR  WILL  B E  USED » APPRGXI MATEH 

THEORETICAL  VALUES  OBTATNFD  IN  PRCORA*  A  . 


VARIOUS  CASES  CP  DYNAMIC  MAKING  A*-1  E  INCLUDED  IN  TH'S 

i^CiS^' — . — AND  4><E  TMF  FOLLOWING  : 


i.-or^A'-'ri"    brakes   u^iriu   a  trnrfrtre   va";aslf 

RESISTOR     . 

.?.-OVNA>i'C    EjRAKTNG    LSING    A     FIXED     gEJjTJilOE 


3. -DYNAMIC    BRAKING    US  (KG    A    MECHANICAL 

AIR     TYPE     BRAKE     . 
-OYNA.'-UC     *R4KiNG     US7^0     A     COMBINATION 


OF    A    VARIABLE    RESISTOR     AND    THF    ^ECHANTCAL 
AIR     BRAKE 
'j. -DYNAMIC — SKAK'TN"     USING     S    C  0  ,J'P  '  'v  AT  V  ON 


OF    AM    AIR     BRAKE     AKC    A    FIXED    RESISTOR 


WHAT     THE    USER    HAS    TO    DO     TS     TO     SELECT 

'NPri;,u;Nr    it   by    tts   mo^ber,   above.... 


A     MtTHCO     BY 


FGR    EXAMPLE 
'/(ILL    BE    AC  CO 
Pr  5!  "iTPK 


ASE=?  ^EAf\S 

MPLTSHtrO     BY     THE 


THAT    DYNAMIC 
USE    CF    £     FIXE 


RRAK  TN 
D 


TN     ADDITION 


LSCR    HAS     TO     SPCCTFY     THE    FGLI 


JVG- 


(A 


)      -     NAX 
CIRCU 


MUW 

T 


1    fi  T  T  C' 

? FS'STG« 
(BJ     -    VALUES 


PERMISSIBLE    CURREN 
<  T  >-I  S     IS    ONLY    LSFP 


r  i    T  Y  F     A  R  M  A  T I 
THE    CALCU- 


f= — TM1!     TIHLOtJ  fcrfft  AL     VALUES     OF     TH>-     V^TM'[[ 


) 

OF 


the  FLFCTkIGAL  CONSTANTS  GF  THE 


'CTOR  Km    GbN':RATOR,KM  A\G  KG  *: "  S  FFCT  i  V  ET  Y  . 
(G)  -  VALUES  OF  THC  MECHANICAL  CONSTANTS  OF  MOTOR 

>MO  GFMFR,ATnlv,.K,MM  AMC  «GG  P  i  S  °E  C.  T  I  V  El  Y — . 


AMQ  Gf 


(D)  -  VALUES  OF  TO 
M\C  GENFRATGR 


:  TNTERNAL  PqS" STANCES  CF  MOTOR 
R,M  ANC  RG 


t*r 


1    VAI  OF  S    FOR — f+tf — TivERTJ.  a    GF    PRUPeL'.  r3 — ; — 
SHAFT     AND    MOTOR     SYSTEW     AND       VALUES     FOP 
fNERTTA    OT    GENERATES     RCTGR     .POWER     TQRn 
A    LTF     I  H  E     I H  E  5 1  S . 


THE 

ME    SYS 


EM 


as   rn    APPEND. 


^TERMINATION    ft    TPEN    LOOP    RESPONSE    OF    THE    SYSTEM 

FT  ho     T'<APZ 
■ T  0  L  H    NPI  OT 

Acr-p    r  a  c  c 

INST    NPLOT=l 


120 


DETERMINATION    OF    PROPELLER    PAR  /i?>i  E  T  tf  R  S 


I L£l£RJLLhAUim    Of    SKLE    EARAMELEJiS 


;CNST    M  =  ?."73C7E  +  C5 


DETERMINATION    OF     POTOR    CHAP- AC  T  EP  1ST  T  CS 


! Pr  tpp^-natthm     PF G£Jl££  ATDR     CJlAB-4£XEBXSXlC-S 

INST    sG=0.2  5  666£»KGG=0.02S6697tRG=C.00C5 

;■  ■■— '  ■'■        ..  1.1,1—     i    ■      .I,  ■ .-   ,.  . — __ _ 

INERTIA    CF    GAS-TUPfUN'E/'oENEfcATCR    SET 


nnn — i'^jm'j. 

I r-THFR    CXN-S1AMXS 


,CNST    P  '  =  1.1  M  5<~  ,P=2. 

.INST    LaSL-1. 


.MTIAL    CCNCITIGNS 


[NCCN    a  C 1  = : 


2  .  0  3  3  3 


0.0,IC4=60  . 


f.P  "VAT'VF 


L-fNTGra  (  IClfUOCT) 


iGT 
TE 
i»T 

P= 


NT 


TT 
Ar'. 


GR 
NT 

■•( 


L(!C3tUJ 

GRl  (  TC4,GTCCT) 


V) 
V  )  *  V 


-E4 


V) 

+  (N*D**2 

■A-J 


=  N 
T= 


ts~s 
P  = 
PR 


TA 
CT 


/8 

>C 


T(S fVP) 

»  (  r  ;>  i?  g  ) 

T  )  *  T  A 
=£1  WL 


(  1  .- 
_LTP 


UDCT=TPR 
V=U 

C^-T/l-iCC  (■>  ,  VP) 


>=c< 


;rl{ 


I*  -  K 


*  *  3  I  *  A 

TC2iMDQT) 


F  N1 -I- fv  o  E  Al 
j-RPJ 


=  K 
=  E 


'J'1" 


TFG 
(  RM 


G=- J  A  *T  F  N *KMM 


iv  = 
AX 


(  7  A  *  *  ?  )  * ! 
-  F  M  *  T  A 


;2|   T;[  »J|-T  fl  ■',  R  M 

ERC'l  =  (  PPB/PMAX  )  *100 
C2  =  (  PPlV  /F"AX  )  *10Q, 


7  MF 


^VNA^aC 


^  (k.i.  !  .0.2166)   ;;i-   rn  run 

LF(CASE.EQ.l)    GO    TO     100 

'Fir^E.EC.')     GO    TO     ?D0 


IF(CASE.EG.3  )    0C    TO    300 
IF (CASE. EQ. 4)     GO    TO    300 
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iFCCASe.EC.S)    GO    TO    3  00 


ma 


JE&s££L. 


[FG=-'+0. 

RB=-0.0093  53*nMf:+Q.017 


GC  to  i2C 
iiO  K«=0.0015 
T7U    UO  !  •  ' -I ' )  h 

GO    TO    60 C 


TO     110 


?C0    SB=0.017 

[Fi"=4Q.  

II-1,  =  40. — — 

C*OT=0.0 


Si    CF    AO    AIR     CLUTCH     BRAKII 

pl=i>.+  rrzrz*?*  (^"-;-m — 


IG    OF  VICE 


TOtr 


P  S  =  5  .  *  T  T  v 

cr.l  =45=5  0.*  (PS-Pf.) 


QVCT=0.0 
IF(OCL  .LT.i 

!&C Til      ?,(-),"? 


.C) 


TC    301 


301    GCL=0.0 

IZl    CCNUNUF 


FWgT=TW  . 

DC!  V='«-550.*FNET 
iFton  .r,T.cn  v] 


G n    TO    3  o  • 


303 


BO'i    CC 


GC    TO    3  04 
CL=OCLM 

TTMIIF 


SF(CASE.EG.4) 
iF( CASE.EQ .5 ) 

■;;r    pq   bno 


GO 
GG 


TC 
TO 


10n 
200 


60U   »SUMQ=OMnT-eP-GFR-QCL+Q.MG+QMGG 

td   ooo 


■m 


5UM( 
GC 

■  c(CASF.qg  .3  .flhO.CASB.LT.B  )     GO    TO    ?U 


=  0. 

TO 


TTT 


CCL 


o 
Hi 


K  = 


t  rM 


777 

F.-6 


.7 


GCl 

IF( 

QCL 

CCN 


L  =  t:> 


50. 
.IT 


=  0. 
T  T  N 


0 
UE 


*  (  F  S I  -  P  C  ) 

.CO     GO    TO    701 


701 

7J2 


TTT 


~U*7 

PSR 

POM 


=TS7 
0. 


PEP 
PEP 

EV  = 
RR  = 


01 
C2  = 
ftnr 

0.0 
0.0 


0.0 
0 .  0 


IAS 

CAS. 


T=-330CC. 
T=/MV» TFV* [AST 


3CC 


GC 
SUV 


TO 
0=0 


ROC 
AST-CP-OFR+QCL 
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vjOO    CCh  T7NUF 


ij  f  ;< 


IT: 


,1'M 


pp. 


P- 


rtL 

1= 

-4-» 


=  60 

B  it  M 


W(N 


,-22500.  ,0.0,  4-225CC.  ) 
?-*P  LJil  I 


/+G 


COCO. 


P  w  /\ 
t  fl 


XI 


P  v. 
A* 


Ay/ 
o.c 


J  fl77THT0 

ore 


LGOOOO. 

1515 


tttht 

CP2 

_k£L 


P/l 

II 


Hi 
C. 

conn. 


+^rc : 

CTRL    FTNTT  'A  =  25.  ,OELT=.l,nELS=.  1. 

;M    .It  FN  ,  PR  B,VO  IT,  OP,  V,  !  A,  PPM  ,  PFKCl  ,QFR,NR  FA  L,  PR  ,  P^«X,  PE  P  C  2  ,  CWC- , 
"  ST,ivAST,Q^'G";,.r,U^.3,S  .Lb, CO, PS  ,OCL  ,03  ,"f  F,C  T  ,  P  C  ,Q^E^  ,'^rr  ,  , 


I  F  G 1  »  W  F  ,  I  A,QE,VP,W,T,A,B,TPR,K,PS( 


THF.  FCLLCWING  STATFMENTS  ARF 

r)i_r)TTJ>iir    PACKAOF     FTif.:    Li  S  L  /  3  r.  0 


APPL>CABLC 
■ST.vui.  AT  <0r>. 


TO 


IF 


l.\\r.\j,\C.F 


AT     M  .  P  .  S     .     T  h  E 
OTHER    TFAN    THI 


STATEMENTS     SHOULD    BE    CONVERTED    7 
;    PLOTTING    PACKAGES    APE    AVAILABLE 


CALL    DRWG(lfl,TIME,VJ 

f.  A  I.  I.    d°wp<  2.1  .i;vr.Nr-F>\L ) 


CAl  l.    DRWG(  2  ,  1  ,T  t.ME,  PPiil.  ) 

CALL    ORWC-l 
ua — awe 


-tw 


-W£-<- 


,2,TTmf,drmi  ) 
.2  .  TI  Vt.  PF.PCl  ) 


CALL    DR WG ( 2  * 4 » T I M E » PERC2 ) 
CALL    DRUG(3t5tTIMF,PWAXl) 


-arr 


iALL    CTrKtrnrrn 

CAL!     QRWG(4»2tTIME,EM) 

CAl  1      npy:0(5.1.Tr^.0r.Ll.  ) 


CALL    EMDRW(NPLOT) 


CP 


TT27T 


TASULAT20N    OF     SPFET    REDUCTION    CPEFFTCTENT     VS. SPEED 


FUrCTTQM     TABMV) 


TTFTF^ 
CATA    VT 


joa: 


l>:  vt (i3»  »ht( nn 

/  0  .  C  ,  ?  .  5  ,  5  .  ,  7 .  5  .  10  .  ,  1  2  . 5  ,  1  5  .  ,  I  7  .  5  , 

r  .  .  q  -■-, .  .  a  i  ,  .  3  2  .  .  Q  ?  3  ,  .  C  4  5  ,  .  0^5 


T/6*C 


-j- 


'.O.t  ?2.5,  2' 
.02  3  ,  .0  2. 


.,21, 


!C.  / 


I  F  (  V , 

BELV=2.' 


30.)     CO    TO    I 


f,-1  n  X  (  V/CELV  )+i — 

SLCPW=(WT(N  +  1)-!nT(N)  )/(VT(N  +  l  )-VT  { l\  )  ) 
TAP^SLfPh*  (V-VT(h  )  K  -J  (  N  ) 


"PTTTF 
T  A8''*=0 . 
■■frFTUN 


0 


ENC 


TABULATION  OF   SHIP'S  PEST STANCE  VS.SPE^P 


Flir-CTTON     TflBR(V) 


DIWFf.S^ON    VTU3)  .P.TTU3)  ^  00  „_    c     „. 

LATA    7  T/ O.C,  2. 5,  5.,  7.  5,  1.0.  ,  12  .  5  ,  1  5  .  ,  17  .  5  ,  20  .  ,  2  2  .5  ,  25  .  ,  tl .  5  ,  sO  .  / 
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fc)ATA    *  TT/ 2*0. 0,7000., 13000., 25000.  ,39000., 570 00. ,00000.  ,10*000. , 

1  i "•coc'o. ,  ^  7:-coo.,  -^snoo . , 2 3 0 c 0 c . /  " 

-  ,n  \/.i,.i  rrcn  en — 

CEI.V=2.5 

n  =  •  f  •■  <  M/  /  c  E 1  v  )  i  ' 


rr 


SLHP^=(f<TT(\  +  l)-RTT(N]]/(VT(K  +  l)-\,T(l\)) 
TAftP  =  «;LCPP*(V-VT(M)  )+RTT(N  ) 


TAnP=20000C. 

Ietupm 

tkit 


TABULATION  OF  THRUST  REDUCTION  COEFFICIENT  ^.""FFC 


FUNCTTGN  TAPT(V) 

■i, j  .-''k.'.ti'im  rrTTrrrrrrrr) 

UATA    VT/O. 0,2.5,5. ,7. 5, 10. ,12. 5, 15. ,17. 5, 20. ,22. 5, 25. ,27. 5,30./ 

r,t\T.\    tt /<:*c.c...ca  .  .vs. .  ■)?.  ..ig  .  .n7c..o??  ..  T^..n?/ 

(  V  /  D  E  L  V  )  + 1 
T.3  0.J     &fi     Tf:     I 


DELV=2 
N=IFT\ 

rR-Vn 


|LCPT=(TT(N+1)-TT(N))/(VT(N+1)-VT(NJ) 
|ABT=SLOPT*(V-VT IN) )+TT(N) 


TAbT=0 

RFTbf<N 


.0 


ENr 


TABULATION    TF    THpiJST    COEFFICIENT    VS. SECOND    vCn"cTFC 
ADVANCE    COEFFICIENT 


FliNCT^      TAPfMS.V/P) 


EI  mpns CON    ST(21) ,CT1(21) ,CT2(21) 


-in 


4- 


.5.1./ 


.1,0.0,.!, •  2  » 


DATA    CT2/-.4,-.15,-.05,.  )5,.1,.15,.2,.26,.36,.4^ 

1.45, .'+2, .39, .37, .32. .13, .50/ 


,  .30, .31, .4, 


Trm4 


A    r.  T  [  /  -  .  q  ;  - .  g  ,  - ;  ?  3 ;  -  .  1  ?  ,  -  .  gg  t  -  .  3  8  ,  -  .  -'t ,  - 
§5,-.34,-.25,-.2,-.13,-.05,.Cl,.10,.21,.45/ 

is=o. 


tr**: 


-^-&- 


-  -j  »    • 


N  =  .^\X(  {  <;+l.  )  /DEL  S)  +1 

IF(VP.LT.C.C)    GO    TO    I 

SI  rr,  Tl,  -  (CTl  (^-t-D-CTI.  (\)  )/(  ST  C\  +  H-S,T(N  ) 


TA. 
GC 


iCT=SLCCTl*(S-ST(N) )+CTl(N) 

Tn   6 


■sloc t  :  -  ( c  t  •»  ( i\  *\ )  -c  '  ?.  <  ■'■:  n-rt-3 

BCT=SLOCT;*(S-ST( N) ) +CT2 ( NJ 

TIJPN 


('■:)  )  /(  ST  (Pvtl)-'ST(N)  ) 


TA 

3F 
TF 


TABULATION    OF    TORQUE    CCEFFIC 

ADVANCE    COEFFICIENT 


VENT     VS.SECONC    vcctmpc 


.funct;dn    tabcq<  s,  vp) 


'''IMhUbi'ON' <T(2ii  ,C'JU21  J  ,CU^(LlL 

DATA    ST/-l.,-.9,-.8,-.7,-.6,-.5,-.4,-.3,- 


■■-U4 


1,0. 


1,  .2, .: 


1 , .  <i>  t  •■<  * 


U- 


CA^A    3Q1/-. 04, -.045, -.05, -.05 5, -.0  6, -.061 
L-.  017, -.05 5, -.0  5, -.0  35, -.0  26, -.01 5, -.003, 


r*f- 


'32/- 


,i3,-.03  5  ,-.  Olrj,  O.J,  .01  ,  .  J.j,  . 
1 .  06  3  ,  .  0  7  ,  .  C  6  2  ,  .  C  6  ,  .  0  5  3  ,  .  0  5  3  ,  .  C  5  ,  .  0  8  / 


,  -  .  G  6 
.005, 

4,.  04 


2,-.C6 

•  J  c  t  •  ^ 

?..0S  , 


4,-.C5 

3?  ,.C7 
rfr?  t  re 


/ 


P I  V P . I T . o . 0 1    on    T 0    2 
-SI  nfYU-lCCI  U  H  i-CQl  ID  )  M  SJ  UM4  )-ST(  M  )  i 


TABCQ=SLOCCl*(S-ST(  N)  )  +-CQ1  (N  ) 

GC    to    6 
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LN^ 


)C.Q2~  (CQ2  fN+l 
>CQ=SUlCg2jMS 

Tnr^ 


-CQ2  dJ)  )/ 
ST(N 1 ) >CQ, 


ST  (N+l» 

(N) 


•  S  T  ( H  )  ) 


TABULATION 


IF  PARAMETERS  OF  TH<=  PRIME  i^PVF« 


FUNCTION  TAEQE(V3,WF) 

DIMENSION  V3T  (5  )  ,hFT  (  10)  ,QET  I  ICO  ,  100) 


)  =  J  .  L 
)=100 

)  =  ?or 


)  =  3  0  C 

)  =40C 

.a 


C. 


C. 


)  -  3 :  i  c , 


)=43C 
)  =530 


)  =  b  jl 

)=73( 
L=a2£U 


7) 


30  C 


8) =10300 
)  -lUrr 


0  )  =  1 2 

,  1  )  =  2 


3  00. 
CCCO. 


■=rr 


=_± 


t  - 
,3 


7  ceo. 
ism 


Jrr-f 


=  5 
=  6 


;500. 

2  C  0  0  . 
e  r.  o  a . 


»  '-5 

.  Q 


=  7 
=  .3 


7  COO. 
f C0G. 


t ;  i"  r="??troo 

2,1)  =1.2  500. 
O     ?  i  —  7  H  0  / 


,  3 


=  2 

-  "3 


7CCC. 
4  5  0  0. 
15C0. 


.7 


=  5 


SCOO. 
4CCC. 


•ere. 

'» C  0  0 . 

gp)^ccoo 


5 


8 
=  1 

=4 


CO, 

3000. 

,-5  0  9  ■ 


k  ceo. 

9000. 


3 '40  CO, 


,  G  J  -  3 
,7)=3 
♦  3)  =4 


Sq00. 
50CQ. 


,10) 


4  ,  1  )  -  5  C 


CCCO. 
5  5  00  J 
r  r  •  /] 


f?.) 

,3) 


cog. 

c.  C  G  G  • 


)  =  \i  :■  e  i 


,5)=2 
.5)=*» 


2CC0. 
6C00. 


cooo. 

4CC0. 


♦  ;  7  1  =  3 

V,,«)=3 

i  ,  *)  -3  8  50  J 


10)  = 
,1)=3 


42000 
COO. 


0000. 
3000, 


,  3  )  =  1 
-4J_=1 


■»5J  =  1 
,  6)  =? 


7T7TT 


4^2-4 


0  000, 


6) 

0)=3 


7C00, 
C5C0, 
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GET 
X3  = 


C3  , 
V3  . 

7T7 
V3 

XLL 


10)  =24000. 
I.  T  .  C  .  C  )  G  C 
GTTTnJTm    G 


3  = 

G 


0.0 

TO 


JLO-L 


100 


G 

B)c  i 

■ r 


TO 
\ti  F  • 

TFT 


100 
LT.?3C0. ) 

rTnrrrr 


GO    TD 

- rru 


T^ 


I 


UF1=W! 


_U_ 


m    ?r; 


23 
TO 
-4r3- 


30. 
200 


20  J      I: 
J: 


FIX(  X3/1000. )+l 

FIX(  I  WF 1-3300.)/ 1000.  )+l 

un  =  ^s-v  jt;  'I  ) 

DW=kFi-WFT(JJ 

Q£j  n-  =  f  hr/ioti.  i* 


OcT(-  +  1  .,)  1-uFTI  T..1  )  ) 


QW= 

CF  = 


(nw/lCOO.  )-MOFT(t  ,j  +  l)-QET(  l 
OET( " , J  J+DELQR+OELQW 


J)  ) 


;no 


fl-P     mtrCTTTTTG    rAT*    APT     HAKT    IF     THE    PLOTT'"!'!    pACH^,e 

cn'>    DSL/360    SIMULATION    LANGUAGE    AT    N    P    S    .     ThEY 
SHQLLLC    EJE    r-''irTF  tf'}    v^    otheh    than   this   PjLCJjriliG   PAjLKAJLE 
TS    AVAILABLE 


'PIC1.  SYSTN     PD     * 

.fit    CF    THF    VFJ_CCTTY    VS.TIf'E 

,HAM  TIM'  g. r? 

,0  6.  -2.  7.  5. 

TT     f.F     IHE     V/\PT/>TTnN     QF     PSHPFM    FR     SPFFP     VS.TTMP 


■  FLAMIN'S  LA 

,C  4.7  -230.  IOC. 

.CT    CP PBLil.  t  pa,vi  ,PEKC1,PF^G2,PMAX 


5. 


FLANT7MS  1A 

0.^4  CO 

ct  cf — i  a  a'vc  ei»  vs  .t;i>*e 


80. 


5. 


1A 


C.C 


.FLA  i\  TINTS 

'.CT    CF    Rb    VsItt^F 

,    FLANT;KrS  1A 

r€ 3^X8 O-UJ- 


150. 


5. 


C.004 


-^ 
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FXE'C    nsi. 

iSL « :rjp' jt  jo 


BJ2XXL&ALI    'Jflr.    !aiLLXI£jN 


LLSJ    UlLQ     ST  VIII  ATTiV 


LANGUAGE  .  TRAMSLATrnN  OF  THE  PROGRAM  WOULD  BE 
HEQUIREH  TF  ANOTHER  SIMULATION  LANGUAGE  TS  USED 


TT 


R  JNRA" 


IHE    ECU  rwT.'j'-.    pphgram    LALLLdLAIES    thp      hppn      i  pop    Pif_££Ciil£ 


nc       Tmf    SYSTEM    UNOFR    DYNAMIC    BRAKING    CONDITIONS. 

THE     TNFRTTA    r]F    MOVING    PARTS    IS    P.EING    USED    EC 
nVNA»IC     RPAKING» 


s^?    has    in    specify    Tmf    EILLLiMIM) 


(  A) 
4-H4- 


-    VAL 
WCTC 


UES    0 

R     AND 


VALUE",    I. 


f    THE 

GENE 

G    tup 


ELECTRICAL    CONSTANTS    CF    TH' 
RATOR»KM     AND    KG    RESPECT! VEl Y. 

■1ECHANTCAL    COM  ST  AN  TS    TF     MOTOR 


(C) 


AN 
-    VAL 


D    GFN 

ues  n 


ERATO 

F     J^z 


-    VAL 
S_bA 


;g\PRATr; 


r?,K^v    AND 
INTERNAL 


UES  F 
T  AN 


OP  TH 
0  MOT 


EM  A  N  C  KG 
E  INERTIA 
Ok  SYSTEM 


KGG 

RES 


RESPECT 

STANCES 


VELY 
n.  F  ?^  0 " 


(D) 


CE  PROPELLER  , 
A\0   VALUES  FCR 


TMT 


INERTIA    OT    GENERATOR    RCTCR    ,POWE 
AS    IN    APPENDIX     A    OF    TEE    THESIS. 


R    TTjRBTNE     SY^tpm 


C Crrrr?-=T! — ;  •  m  ;  t  !  fl  l_  y  5  E  S — l"HE     "HASE-ST'JP     CO'^maMJ 

CRSTP=1     ,     PRGHI RITES    THE    APPLICATION    CF    THE    EMERGENCY 
PROP.EDURES    ANi.     TEE     S  E 1 P    OPERATFS        AT     COMMIT TCN 

REGULATED    BY     THE     INITIAL    CCNDITTPAS        AND 

STEADY    STATE    DYNAMICS* 

T6G    TRAPZ 

GER     NPLQT 

T.GEk    CPSTP 

VST    NPL'jT=l 

r,ETFPMTNflTir)N    OF     PRPPE'LER     PARAMETERS 


NST    1=33171.  ,C=13. 


-•  t 


DETERMINATION  OF  SHIP  PARAMETERS 


viM  M  =  2.'i'3^7F  +  C5 

.;ETEPMT,\ATIO\    CF    i'OTQR,    characteristics 


VST    KM=O.l63  26,K»>v=0.4  59l53,RM  =  G.CC05 


DETERMINATION    OF       GENERATOR    CHARACTERISTICS 


"Si    K(!  =  o.?366A/i,KGr,=  0.03«(8>6'5^R<;  =  C.'JuCb 

— T^PRT7A    OF    GAS-TUR^JMr/JC:IERATCR.    SET 


*.ST 

IG 

=3916. 

OTHER    CCN 

STANTS 

■  '■j  SI 

PI 

n 

",  T  "  A  L    C  0  |\  C 

=  2. 
VT 1 0  N  s 

:CN    iCl  =  30  .  ,  I C 2= 3. 8 3 33, I C 3=0.0  ,IC4=6G. 
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ST 


,TP  =  1 


IM2  ~'" 

U  =  i  Mir  i.  (  ;n.  ,uoct) 

X=- NTGRL ( TC3,U) 

DT^iVlii  (  V) 


A8W( VJ 

(  ]  .  -  W  J  *  v/ 


-1  =  I  .vO1  (  v ) 

p  =  nr,r  t  (  A  ) 


CT=TA8CT ( S  ,\P) 


TP 

t  r  i 


(1  .-T  )*TA 
=  ( TP-RT  )/M 

ra  rrp 


-urr 
r.c  =  T/tf'.r.fj(  Sa^EJ 


£P=CG*P*(n'S'*2)*A 

Bgt=:ntgrl ( TC4,nTonT) 


|FM=40. 

IFG1=(  EM-32.  )/(KG*NGT) 


1  [i- itf  b  ;<,  w ---  ,'  \-  V  ?  ,'v  P  1^  A  L — 
EG=KG*" FG1*NGT 
[A=j  F'".-Fv  )  /  LE£±fL! 


W-T&- 


TF (CRSTP, 

; f ( r:r  p,tp. 


ECO.)     GC    TO    110 
n;.l.AM:.i'i.tT.o.g  ttns ) 


GC    TC    10  J 


Lwf  =  : 


.0  0 


GO 


GC    TC     ?00 
WF=3? J0.+1CCC. 
1 F  (  i.i  r .  r. 


sR. 


FC.if  ,r,E.l"?00.  )    Ci 


MP(1.4) 

c. 


t-U^- 


CC    TO    12C 

feF-1230G. 


<i(j    CENT  'M1IF 

CC    CC'\'TT\jiip 

g  y  =  k  M  ■*  ~  T  A  *  T  F  W 


cc-FM  =  Knc--  t  <i*tfc;i 

GF  =  TARQC(N3  ,ViF  ) 


GTCCT  =  ntp/  (? 
N3=.\'GT*6  0. 


*PI 


^^C-OV.-CP-CPR 

EFR  =  FCN!SV»(N',-22500.,0.0,+2250C.  ) 
lDCT=SUMC/ (?.*P1*I  ) 


i\  f*  b  A  [  =  t-.0.*h 


TRL    FTM 

frT    .1  ,5,V,QP 


.,CELT' 


C. 001, nELS=J. 0125 

T.CTf FT  tC'M,  latf,  ;  CC^QC-rM,  Qf  ,u|f  ;  A,T,TPR  , 


GiF.!lM^BlTA1iTP,N3iEG,SUI;1Q,0FR,NDCT,NREALtXtIFG2. 


:FG,«CtPF,EF 


THE     FOLLCVTMG     STATEMENTS     ARF     APPLTC/-6L£     TO     THF 


PLOTTING  PACKAGE  FOP  OSL/3fcO  SIMULATTC\j  LANGUAGE 
AT  M.P.S  .  THF  STATFMCfvTS  SHOULD  BE  CONVERTED  IF 
OT'IFR     TI-ftN    THIS     PLOTTING    PACKAGES    A31  E    AVATLAuLG     . 


CALL    ORV*G<  l,l,TiMb,V) 
CALL    OPKG( ?  ,1  ,TT ME, NREAL ) 

cat  l  ^py-.Gt  ;  ,i,tivf.  TFr,i  t 


CA{  L 
CALL 


DRWG( 


FG) 


DRWG(4»1  ,TIME,N3) 
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LAll    0PWG(5,1,TIME,RC) 

C  d  L  L    !'.  P  I* r-  (  5  ,  2  ,  T  ;  m  E  ,  g  H 


M  1^  A  l 


rNI)H'»  [A£l   nT  1 


T&AN 


TAi4ULAT(PN    OH     M'FEP     PEnuCTTCN    COEFFICIENT     ^S.^FEn 


PPWtTfOT     fN.'MVI 

BIMFNST CN    VT  (13  )  , WT ( 13 ) 

1 1  ,a  I A    \/  T  /  ( 2 .  G  ■  2  •  5  t  5  •  t  T  •  '5  .  I ' )  .  .  1. 1  •  5  .  1 J2  .  j  1 7  .  5  ,  2  0  . ,  ?.  ?  .  5  ,  2  5  .  ,  2  7  .  5  ,  3  0  .  / 


DATA  WT 
TF(V."T 
4-^ 


/4*C.C,.OQ5,.01,.0  2,.C!23,.l.45,.C45,.G-:!e,.02,.OC'</ 
,3C.)    GO     TC    1 
s 


N=7F iX(V/CFLV)+l 

51CPW=  (  *T(N+] .  )-WT(Nl  )/(  VT(  N+l  J-VT  (M  J 

T  A  K  W  -  3  L r  PN  -'  (  V - V  T  < \\  )  hFWTI  M  ) 

RETURN 

TAB>;-0.0 


RETURN 

ENn 


TA«LLATft:r'    91s SUP'S    RESISTANCE     VS.SPEFC 


FUNCTVQ 
ll'WF.N.S? 

OAT  A    VT/Q.C  ,?.5 


N     TAPR(V) 

ON    VT(13)  ,RTTU3) 

5  .  .  7  .  5  .  1 0  .  ,  I  2 


-3-* 


15.  ,17, 


-23- 


'J^J- 


DATA    RTT /2*C.O» 7  000. 1 13000 •  ,  25000  .  ,39000. , 5  7000., 90000 
1136000. ,  173000. ,225000. , 280000.  / 


,10300  0., 


J!P(  V.n,T  .30  .  )     GO     TO    1 

BEl_V=2.5 

IS.=  T  F V  X (V/HFLV)  +  1 


SLL  PI-  =  I 

|abr=sl 

SFTIiRf 


k  T  r  (  n  + 1 )  -  p  t  t  { r; )  )  /  ( V  r  (  n  + 1 )  -  v  r  (  m  ) 

DPR* (V-VT (N) ) +RTT(N  ) 


TAP.R=28 
PETURM 
-trtrf) 


00  0 


TABULATE1    OF     THRUST    REDUCTION    COEFFICIENT     VS.SPEEC 


FUNCTION    TABT(V)  

Btt*gtn*Q1^    VT(13)  ,TT(  13) -~-~ ~    "     ~ 71     "     ~     . 

DATA,    vt/o.  0,2.5,5.  ,7.5.10.  ,12.  5,  15.,  17.  5,  20.,  2  2.  b,  2  5.  ,2  7.5,30./ 

DATA    TT/5*C.C,.C1,.Q5..Q7,.C3,.0  75,.072,.06,.02/ 

UbLV-2.  5  — 

N  =  IFTX  (V/0EI.V)  +  1 

1  F(V.r»T.3Q  ■  )    Gf,    TH     1 


SIGPT=(TT(N  +  1)-TT(M)/(VT(N+1)-VT(N)) 
fABT=SLGPT* (V-VT (N) ) fTT(N) 

"HrTLiSM — 

TABT=0.0 

RETIJKN  


"ETTJ- 


TABULATTON    OF    THRUST    COEFFICIENT    VS.SfcCONC    PCDIFIEI 

APVAi\CF    COFFFIC  I  EM  
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-\  iiM.  i !  rm — '  a h l  1 1  s ;  v p  i 

CI  MENS  [ON    ST(21) ,CT1(21) ,CT2{21 

nf,T".  _ST/-1  .  .-.r,  .-.  9. -.7. -.;■,.-  .' 

1.4,  •  b ,  .61  .7  »  •  8 ,  •  9 »  I  •  / 

LATA    CT2/-.4,-.15,-.05,.05,.l,.15,.2».26,.36,.44,.30,.al,.4, 

4  ■  *  5  ,  ■  4  3  1  1  3  g  ,  .  3  7  ,  .  ?■  3  ,  .  3  3  ,  .  b  0  / 


~  *  ■  ' — ■/«~.?i~.1  .  ( ■  .  0  .  .1  .  . 


*Li_ 


DATA    CTl/-.4,-.2,-.28,-.3ci,-.35,-.38,-.4,-.41,-.3  6,-.2c,-.?5, 
L-.3  5»-.34,-.25,-.2,-.13,-.05,.01,  .10, .21, .45/ 


■I.F.I  <,=  3.1 
5=IFIX  (  (  S-t-1.  )  /OLLS)  +  1 

Ff7!'.i.T.c.oi    r,r    to   ? 


Sl.nCTl=  (CT  KN  +  1J-CTH  Rl  )  /(  ST(f\+l  FTtTnI  j 
TABuT=SLCCTl*(S-ST(M)  )+CTl(N) 

-T,C    TQ    6 

|L0CT2  =  (CT2(N+1J-CT2(IM)  )/( ST(N+1 J-ST(N) ) 

TAeCT  =  SLCCT2*(S-ST(N)  1+CT2  (N) 

■"<f-Tij;;,"j — 

END 


TABULATION    OF     TORQUE    COEFFICIENT    VS.SECONC    VCCIC*EC 
ADVANCE    COEFFICIENT 


Ftfl  iTtffl    TAPCQ(  S,VP  ) 

DIMENSION       ST (2 1 ) , CQl ( 21 ) , CQ2 ( 21  ) 
OAT  A    S^/-l  ♦  ,-.? ,-.0, -.7,-. & ,-.5 ,-, 


4t~»3t-.?i-.  LtO.  t  .1 


1.4,  .5,  .6,  .7, .8,. 9, I./ 
DATA    CQl  /-  .04, -.045, -.0  5  ,-.055,-.0  6,-.0  61,-.C62,-.C=4,-.C5>;,-.04  , 


4— 


■  ". 055, -.0  5,-. 335. -.026. -.015, -.00  b. 


005  ,  .02 


C7/ 


OAT  A    CQ2/-. 08, -.03 5, -.01 5, 0.0,. 01,. 02, .034, .042,-06, . 07  ,  .045  ,  .0  5  , 
ti. 068 ,  .  07  ,  .  C6  2  ,  .  0 6  ,  .  05  3  ,  .  05  8  ,  .  C  5  ,  .  0 8/ 


OELn>=0.l 

N  =  TF?X(  (  S+l.J/DELS)  *1 

::f(vp.lt.o.o   go  lo   ? 


gUJCQT=TCf]  [TN+I  r-COl  (N)  )  /(  ST  IN+  l)-ST(N).) 
TA6CU=SL0CC1* IS-ST(M ) ) +CQ1 (N  ) 
GC    TG    6 


SLGC02 

TAP.r.O  = 
riFTl'PM 

END 


ICO? (M+l) -CQ2JN) ) /( ST(N+1)-ST(M) ) 
SL0CC2* (S-ST (M) J+CQ2 IN) 


TABULATION    OF    PARAMETERS    OF    THF    PRIME    MOVER 


FUNCTION     TABQF{V3,WFJ 

DIMENSION    V3T I 5) ,WFT (10) ,QET(  100,100) 


V3T (D-0.0 


V?T( 
V3T( 


2)=1 

3)=2 


J  •,  L  . 

30  0. 


~7TT 
&3T( 


4T^J 
51=4 


WFT(l  ) 


on- 
occ 


UFT( 
*FT( 


')  =  43CC 

3J=s: 


3CC. 
fTK)  -63CC. 


Tf 

WCT( 
'lFT( 


5)=7 
6J=8 


30  C 
3CC 


TFT 
%FT( 
ILL 


77^ 
8J  =  1 


30C. 
03C0. 


i»FT(Q)  -1 11  GO. 


lr*F"T  ( 
QET( 


10)  = 
1,1) 


ti- 


•f- 

3) 

4) 


123CO. 
2CCC0, 


28500, 

3 7 CCO. 

■45COO. 


OFT( 
ccT, 


XFTT 
GET 


YTTT 
1,6) 


5350TJ. 

6?CC0, 


iiFTf  ],  ,7)-7CC0n 


QET 
«JET( 


8) 
9) 


^77CC0. 
86500. 
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I '  f-T 


Jl£I 


wdT 

QFT 


10)  = 

=  1 


S5C00. 
2500. 


=*«■''.'     I.'.'. 

=?7CCO. 
=3Ji5£& 


1  5CC. 

fooo. 
i\  c  c  s . 


OFT 
OFT 


,r> 


r?  Q  0  0 . 
A  C  C  0  , 


■ c^-r 

qft 
oft 


)-gqcuo, 

=8000. 


_l2 


OFT 
OFT 


i?4 


PRCO. 
AC  CO. 

s  Q  q  o  ■ 


OFT 
GFT 


i »  b 
,7 


4CCC. 
9500. 


T7TT 

QFT 

Jitl 


gygt-ffr; 


1= 


C  0  0  Q . 
55000. 


OFT 


t-r 


ccc. 
rcc. 
2  c  g  8 


OFT 
OFT 


QET 

OFT 


,5 


."7*00. 
!20C0. 


^Tf 


CCOO. 

4000. 


,  7 

O 

1 


^TTL. 
43000 

oeo, — 


0  I:  T 


q 

;  i 


)  = 


OFT 
OFT 


,?. 


=  t 


c  C  C  • 
CCOO. 


CFT 
GET 


i  i  ?  c  e 


70^0. 
CCCO, 


QET 

CFT 

-0£X 


77 


QET  (5 
|.F(V3 


1  F  (  V  S 

X?  =  V- 
|C   T0 


tlC 
.IT 
TS-f- 


2 


4CC0, 
7  COO. 


)  = 


24000. 

.C)    r-C    TO    1 

on  to 


e-e 


i- 


100 


X~  =  0.'J 

GC    TO    10  0 
X  ":  -  i.  0  3  0 . 


GO    TO 
I.FUF 


-rF-f- 
|F1=W1 
GC    TO 


100 
LT.3:CC.)     GC 

GT.izrcc.)   Gl 


TO    ' 


200 


GC    TO 

c 


201 


Ucl-12  300 


;cq 


i:=IFT. 

j='  ft; 


(X 

( ( 


3/ 

WF 


QP-X2' 
DFLCF 


-V3T(T ) 


I.CC0.  )+l 

1-3300. ) /1000. )+l 


Y5- 
-W 
(D 


FT{  J) 

P/1000.)*(QFT{ 1+1 , J  )-UET  (  I  ,  J)  ) 
U,J  +  iJ-ULI(itJ)J 


fABGF-- 
BETIJFI 


TTTTT7TTJTTTJTT^TTJF 
QET(  I,J)  +  CELOR+DEL 


}W 


END 


TUT     rOLLCWTNG    [UTA    APr     PAj-'T    OT     TIC    PLOTTT'-'G    PACKET 
FOR     DSL/360 
SFOUIF    eE    W 


S  AVAL LAHLb 


SIMUI  ATION  LANGUAGE  AT  N  P  S  .  THEY 
JCTFIED  TF  CTFER  THAN  THIS  PLCTTTMG  PACKAGE 


I!  n  ,  cy  c:  r,y     pp 


CF     THF 
■LANTJM'i  S 


VELCC 

1A 


TY    VS. TINE 
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b .  -  2 .  7 .  5 . 

'  T    Li-    Tt-F     VAP;/>T^rN    HF    PRi.-1°E'.Lf:K    SPFEC    VS. TIME 

nhi  i  .1^1  s rn — 

C                      '+.1                     -220.                      ICC.                   5. 
T   (f      \>q  -,1  .P9M  .pppr  i  .P^r?  .pviax 


E.  ANT  INT S 

0 . 1 A 


1A 


0.0 


f-f— fc-f ta     /\^r     FM     VS*TIMf 


80. 


MAM-NTS  IA 

0  I  .  0  .  G 

r  I     L  H     ^  K     ^.TT,"F 

FLANT'tMTS                IA 
L 0  .  3  Q 0.0 


150. 


U.Q04 


J2^_ 
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//  Extc  nsi. 

^/DSL.  INPUT  '3D 


>_  d ' !  ^  l  J  A  G  E  .  TRANSLATION  OF  THE  PROGRX^WOULD  BE 
ReOUIREC  TF  ANOTHER  SIMULATION  LANGUAGE  75  USFD 


rrrtt pro spa* — c — — — — 

* ThF     FPl.LrwTNG    PROGRAM    CALCULATES    THE    CI  1  S  E  n,    LQQP    REjPCNSE 

»  vkHcM     THE    INERTIA    OF    MOVIMG    PARTTS    OF    THE    GENERATOR 

f A  MO    TH-     PC^ER    TUtftUNE    A-fti    USED    FCk    DYNAMIC — &R-A-K4^1  g  . 

>    : 

r " 

i Tri''_     USrP    HAS     TO    SPFCTFY     THE     FOLLOW, NG     : 

r (A)    -   VALUES   oTF    TEE    FUFcTRTOt   constants    rr.    th- 

:<  MTJTOR     AND    GpOF.<  ATUR  t  KM    ai'nC    KG    p  ESPFCT  *  V  E!  Y  . 

j (  D     -     VALUES    OF     TF'E     MECHANICAL    C  0  \  ST  A  fvi  T  S     OF     MOTOR 

>  AND    GENFRATOn  ,KMM     ANC    KGG    R  LSPbCT  T  VFL  Y     . 

>  (C)     -    VALUES    OF    THE     INTERNAL    RESISTANCES     OF    MOTO& 

> ANC    Gr-^F^ATO1'.  t     *fi    A,\C    ?G 

>  (0)  -  VALUES  FOR  TMC  INERTIA  OF  pwCPELLcR  , 

> SHAFT  AND  MOTOR  SYSTEM  AND   VALUES  FQc  THF 

5  "  T'lHPTlA  m  C'-NFRA  tTR  RLTfjR  f  PUw'rR  TLk  d  I  N  =  HTI  u '" 

*  AS  IN  APPENDIX  A  QF  THE  THESIS. 

_       CPSTP  =  c  ,  TN'  T:  AL'f  SCS  THE  CRASH-STOP  COMMAS 

*  CRSTP  =  l  »  PRQHTRITPS  Ti-c  APPI  ICATION  Oc  THC  EMERGENCY 
! PPIiCt'ntJ°rS  AND  THF — g  I-  [  P  OPERATES AT  COM™  IT  5  Of'  S 

*  REGULATED  3Y  THF  INITIAL  CGNDI T  ttas   AND 

I STEADY  STATE  DYNAMICS. 

c 

I _ __ 

(MTEG    TRAPZ 
[NTGEK    \PLHT 

NT^K    Cr-'STP 

7JNST    NPLOT^l 

3 H.-  n-PK  r  NATION!    CF     Pk'OPELLhk    PARAVbTERS 

'.'INST    I^im,  ■  ,0-1^.^4 — 

s  DETERMINATION    OF    SHIP    PARAMETERS 

► . 

;ONST    M=2.7397E+C5 

- Ot:  I  r:PM  "NA  I  aON     OF     MOTtlR     LHA^ALIbKJSV-'CS 

DE  TERM  'NATION    OF       GENERATOR    C  H  *R  ACTEP  T  ST  *C  S 

JONST    KG=3.2  56  66  6  .KGG-3.0?  5  6  6°  7  ,  RG=C  .00  C5 

r ;MirRifi     IF    GAS-IUkhiNE/GENbkA I L*    ST1 

T,\<;t     ir-i^K, 

> 

'  OTHER    CONSTANTS 

* ■ ■ 

:PNST    PI  =  ?.14J.5^,P  =  2. 

'  i  N  i  i  I A  L    C  U  N  C  .1  r 1  l  ]  N  s 

'i 

Lmt.n    !n^i..Tr7^..qrv.l  r^=n.n  .  To.4^n. . 

XNST       CKSTP-=l 


133 


CRAr     mim  i)  ,nbw(?  j  ,cu  l  ),f  2(  i  ) 

riLF      M'M1  U  =  H.7feF  +  o<nnE-l(l-3)=l.  ,  30  0.  ,19.  76E+04,  C  I  ( 1 )  =10  .  ,C2  ( I  )-40 


j;  '■' '/  ■-> "'  f  y  " 


tMNTGRt  (  TCI  ,UHCT1 

X=fNTGRL ( TC?,U) 

(  T  -  T  A  >\u  (  y  ) 


b=TAB  'J(  V  ) 
VP=(i.-W)*V 

T  =  T  U'  T  t  V  ) 

A*VP#*2+ ( k *n )#*2 

ir.snr  t  (  a  j 

S=n*o7F 
CT=TAar.T(  s  ,vp) 

T  fl  -  g  T  ■  t  P  *  (  C  a  *  2  )  *  A 


TP=d..-T  )*TA 
7PR=(  Tp-RT  )  /M 
U  '•.  r?  T  =  T  P  p. 

v=u 


CP=CQ*P* ( C**2)*A 

NGT  =  INTGRL  (  TC4 ,  GTnf'T  ) 

.-j-t r.iT^nt  (  U,?  ,NDQT  ) 


'FP'^40. 

iFG=(EM-33  .  )/(KG*NGT ) 

t:i:  =  H:MSif'<  i  Pfi — ,-h jp. , u.c, +400 . ) 

[T7JT=TF7  k '" — — 

gf=KM*I  F  M*|\R  E  AL 

"iC-Kr.-io"  Fr.y^'iT 


t/k=  (EG- EM  )  /  (RG+RM  ) 


'NAM  to 


*F(CRSTP. EC. 1. AND. N.LT. 0.21661    GG    TO    10G 


WF=3300. 

GG    TO    ,?0  0 
100    '*,?--!  j'J 0.+  1  CCfl  .^ioiMH  (1.4) 


IF(WC. GE. 12300. )    GG    Tn    no 
GC    TO    120 


110    «!;=1?^)0. 
120    CGMTNUF 

2J0    (;CKr 


QN=K MM#I a*TFM 
QGEN  =  KGG*T A*  TFQ1 

'.vr-TAOQL  H  3  fW  ) 

CiF=QE-QGEN 
iTCOT=OTF/ (2.*PI*IG> 


J  =  N'I  *6U. 
SUPQ-Q«-G-P-GFR 

wFP=n:,'MS'i(  l  .-2-5r,0.  ,f!.n,  +  22Sn,1.  ) 


NDCT=SUMQ7 (2.*P 
NftEAL=60.*N 


;'PLf-  _ 

KiUL    t-!  N!  <  M  =  iib  .  ,!)tL  I  =J.  .)■').».  TUEUS^UTTTZ^ 

INT    . 1,S,V,QP , TFG1,NGT,CT,RT, QM,IA,Nf CO»WfQGEN,Q£,WF,A,T,TPR,.. . 

n-P.-  4tr  t  ta  ,tp,m-i  ,  cpt  <;iimo  ,  i  j  E  S  .  U  Q  r  I .  MP  E  Li  .  x  .  T  E  G  2 .  T  E  G .  P  r  .  P  c  .  r  ! 

— TUG     re  LC'-iT^p    ^TATTHnviTS     »ft€     frP-P-fc  I  C  A  i'L  E    TO    T'ir 

PLTTT-N^    PACKAGE    FTR    DSL/360    SIMULATION    LANGUAGE 

it    M.P.S    .    THE    STATEMENTS    SHOULD    BE    CONVERTEO    TF 

T-PTFT — IFAN     I  H  i  S     PL-'llliNG    MAGKAOES    AkE    AVAILAHlE     . 


GAl  L    DRWGI  1  , ltTIMEtV) 
GALL    DRWG< 2 , 1 ,T I ME, NREAL J 
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IM  L    0"WG(  3,l,T-v<F,  TFGl  J 
EAL  l.    DRWG(2»2|TIME,IFG) 


r  J I   l 


H':t    *   ,   L,    I    l**h,U 


"■»  »  I  t  I  I  "":  ,  l\  3  ) 
M_».    ORVlGt  5,1, TIME, RC) 
r,'t  i     Q£Jci£  (  s  .?.n  MF.m 


-"    (,  ALL    ilM'T^  (rjPLOT  ) 


T  r.  A  ,"v 


TArtULAT'nn;    np     SPEEfi     H^DUCTfQM    CO  FFlv  C  !  i2AT     VS.^PrEO 


fuict'On    t^bU(V)  ~ 

BIMENSION    \/T  (  13)  ,*T(  13) 

BAJ-fl    VT/Q.C  ,->.5  ,  5.  .7.5  ,1,1.  ,12.5  ,  1.  S.  ,17.5,20.,  ??.■; L,  ?5.  ,27.5,  20./ 
U*TA    dTM#C.C, .005,. 01 ,.02 ,.033,. C45, . 045 , .03  8 ,,02, .004/ 
[P(V.GT.30.)    GO    to     I 

P  E  L  V  =  ?  .  5 — — 

N=ICIX( V/DCLVJ+1 

s  i.  r  g  b  =  !  a  T  ! N  *  1 !  -  j T « '" !  !  L !  v r  <  N  + 1 )  -  v  t  (  n  ) ) 


T  5 H  Vi  =  S  L n  PU  M  V - V  T  (  \  )  )  f  ■.* T  (  M  ) 

RETURN 

TA^-Q.Q 


RETURN 
END 


TABUL 

AT  i UN    n|- 

SHIP'S 

RtSl STANC 

E    VS. 

SPEfcC 

FliNCT 

■  0  N    T  A  8  R  ( V  ) 
SIGN    V  T  U  3  »  ,  R 

T  T  ( 1 3  ) 

• 

U  4  1  A     i  w  ij  •  l  i  /  .  j  i  j  .  i  i  •  .'  i  l  ;  .  i  k  •  ;  t  w  • 
DATA    RTT/2*C.  0,7000.,  13000.,  25000., 
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